DTP/00/54 
MPI-PhT/2000-35 
TPJU - 9/2000 

",3" ■ August 10, 2000 

O' 

o 

(N 

_ ' ANALYTICAL QCD AND MULTIPARTICLE PRODUCTION "- 

(D 

^ ; VALERY A. KHOZE 

>0 , Department of Physics, University of Durham, 

CN . Durham DHl 3LE, UK 



WOLFGANG OCHS 

Max Planck Institut fur Physik, Werner- Heisenberg-Institut, Fohringer Ring 6, 

D-80805 Munich, Germany 



JACEK WOSIEK 



> 

00 

(N 

(<-^ ' Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, 

(^ , 30-059 Cracow, Poland 

o 

flj We review the perturbative approach to multiparticle production in hard coUision 

(^ , processes. It is investigated to what extent parton level analytical calculations at 

low momentum cut-off can reproduce experimental data on the hadronic final state. 

Systematic results are available for various observables with the next-to-leading 

1 J - logarithmic accuracy (the so-called modified leading logarithmic approximation - 

1^^ ' MLLA) . We introduce the analytical formalism and then discuss recent applications 

jH , concerning multiplicities, inclusive spectra, correlations and angular flows in multi- 

^ - jet events. In various cases the perturbative picture is surprisingly successful, even 

for very soft particle production. 



''To be published in "Handbook of QCD" (loffc Fcstshrift), cd. M. A. Shifman (World 
Scientific) . 



Contents 



1 Introduction 



2 Perturbative Parton Cascades and Jets 



2.1 Jet Definitions and Multi-Jet Structure 



2.2 QCD Cascade and Evolution Equation! 



5 

6 

9 

2.3 Evolution Equations for Jet Observable^ 12 



2.4 Parton Hadron Duality Approaches 17 



2.5 Space-Time Picture of Jet Evolution 19 



3 Multiplicities 



20 

3.1 Asymptotic Behavioi 21 

3.2 Full Solution in DLA 22 

3.3 Modified Leading Logarithmic Approximation (MLLA)| .... 23 

3.4 Numerical Solutions 24 

3.5 Experimental Results on Quark Jets 25 

3.6 Test of Jet Universality 30 

3.7 Comparison of Quark and Gluon Jets 30 



4 Momentum Spectra of Particles Inside Jets 



4.1 Evolution Equation for ^-Spectra and Approximate Solutions 



4.2 Moments 



4.3 Predictions for the ^-Spectra 



4.4 Comparison with Experimental Data 



5 Multiparton Correlations Inside Jets 

5.1 The General Structure of n-Parton Angular Correlation^ 



32 

33 

35 
37 

41 

44 

44 



5.2 Distribution over the Relative Angle - Comparison with the Ex- 



periment 46 



5.3 Azimuthal Correlations to the Next-to-Leading Order 48 



5.4 The Critical Angle and the Local Parent Child Correspondencq 51 



6 On Inclusive Properties of Heavy Quark Jets 



53 

6.1 QCD Bremsstrahlung Accompanying Heavy Quark Production 54 



6.2 Specific Features of Events Containing Heavy Quarks 56 



7 Color-Related Phenomena in Multi-Jet Events| 



65 

7.1 Radiophysics of Hadronic Flow^ 65 



7.2 String/Drag Effect in qqg Events 67 

7.3 Muhiphcity Flows in Three- Jet Events] 71 



7.4 Topology Dependence of Three- Jct-Evcnt Multiplicitj; 74 



7.5 Experimental Studies of Three- Jet Events 79 



7.6 Collective Phenomena in High-p^ Hadronic Reactions 83 



8 Conclusions 88 



1 Introduction 

A characteristic feature of high energy collisions is the production of many 
hadrons. At existing high energy colliders with incoming electrons or protons 
the mean number of produced hadrons ranges from 30 to 50, at colliders with 
incoming nuclei it is even above 1000. For a long time multiparticle dynam- 
ics has been studied within the framework of phenomenological theories and 
models basing on the quark structure of the hadrons involved in the consid- 
ered process and the analytical structure of the scattering amplitudes, see for 
instance Ref. 1. 

Within the parton modelnO each constituent of the hadron carries part 
of the total hadron momentum and could scatter on other partons or interact 
with leptons and fragment finally at larger distances into the final state hadrons 
under the action of the confinement forces. In the multiparticle final states in 
the "hard" collisions which involve large momentum transfers typically several 
jets of collimated hadrons appear. They can be related to the partons emerging 
from the primary hard interaction. The existence of spectacular jets of hadrons 
- the footprints of partons - is among the most striking phenomena in high 
energy physics. Jets fronuarimary quarks were discovered at the e^e~ collider 
SPEAR at SLAG in 1975Q with the angular distribution as expected from the 
production of spin i quarks. 

With the advent of Quantum Chromodynamic£l (QCD) a quantitative 
treatment of many phenomena involving hadrons became accessible. The force 
between the quarks is mediated by gluons which interact also among them- 
selves. They may be produced in a hard collision process and evolve into a jet 
of final state hadrons. The gluon jets have been observed in e+e^ collisions 
by the four PETRA experiments at DESYU in association with the two quark 
jets. 

These discoveries have opened up a large field of jet physics with detailed 
investigations at the e+e~ colliders LEP and SLC, the pp colliders at CERN 
and Fermilab and at the ep collider HERA. It will certainly remain one of the 
main topics for studies at the colliders of the future. In the quantitative treat- 
ment of jet production one starts from a precise definition of a jet which refers 
to a resolution parameter in a jet-finding and jet-defining algorithm. In this 



way the multiparticle final state is reduced to a final state with only a few jets 
characterized by their 4- momenta. In the theoretical analysis of jet production 
one assumes that the jet of hadrons can be represented by a jet of very few 
partons (one or two, typically) and their production properties are analyzed in 
perturbative QCD. Because of the asymptotic freedom of QCDlil the couphng 
constant becomes sufficiently small in hard processes so that reasonable results 
can be obtained from perturbation theory. 

Today many phenomena in jet physics are quantitatively described by 
perturbation theory in terms of only one basic free parameter, the coupling 
constant a^ at a given energy scale; in addition some structural characteristics 
of hadrons at a given scale are needed as well. For many quantities results in the 
next-to-leading order have been obtained. The accuracy of such calculations 
can be illustrated by the error in the determination of as in jet physics which 
is less than 3%.H The accuracy in these studies is limited typically not by the 
experimental statistics but by theoretical errors from the scale uncertainties 
and, in particular, from the uncertainty in the transition from the QCD partons 
to the observed hadrons. This latter difficulty makes it especially desirable to 
study the transition from partons to hadrons in greater detail. 

In the popular^ QCD-basciLMontc Carlo models for multiparticle pro- 
duction (HERWldi, JETSEtI3 or ARIADNeIO) based on the Lund string 
modeli^ or cluster fragmentation the perturbative evolution is terminated at 
some low scale Qo (typically Qo ^ 1 GeV) for a dynamical variable (such as rel- 
ative transverse momentum or parton virtuality). Then the non-perturbative 
processes take over and the transition into hadrons is described by phenomeno- 
logical models. In these models various parameters have to be fitted by the 
data which makes it sometimes difficult to directly trace the connection of a 
fitted observable to the underlying QCD dynamics. Let us emphasize that all 
existing phenomenological models are of a probabilistic and iterative nature 
when the fragmentation process is described in terms of simple underlying 
branchings. Their success in representing the data relies on the fact that in 
certain approximations it is possible to absorb the quantum-mechanical inter- 
ferences into the probabilistic schemes. 

Alternatively, one can try to compare the perturbative results directly to 
the experimental data without a complete hadronization model. In this way 
the number of non-perturbative parameters is minimized (in addition to the 
QCD coupling a^ one introduces a non-perturbative cut-off Qq in the partonic 
cascade). Most importantly, an analytical treatment becomes feasible which 
reveals the details of the QCD dynamics. Also non-probabilistic phenomena 
become accessible. In recent years two kinds of applications of this general 
idea have been developed. 



An application of analytical methods is suggested in particular for observ- 
ables which are "infrared safe" , i.e. do not change if a parton splits into two 
collinear partons or emits a soft parton. Then the dependence on the non- 
perturbative cut-off Qq is suppressed. In the last years various applications 
to shape variables such as thrust and others have been carried out. The per- 
turbative results for the mean values obtain characteristic power corrections 
0{A''/Q'') with predicted power k. Together with the 2-loop calculation for the 
lowest orders in perturbation theory this approach is used to obtain accurate 
results for the couphng asJij 

Another analytical approach concerns particle densities and correlations 
which are infrared sensitive and depend explicitly on the non-perturbative 
cut-off Qo- The concept of "Local Parton Hadron Duality" (LPHdHJ) has 
been formulated originally for inclusive spectra and states that for small Qo 
of order of a few hundred MeV the parton distribution calculated in pertur- 
bative QCD gives already a good description of the observables. Indeed, the 
general features of the inclusive particle distributions in jets produced in e~^e~ 
annihilation, deep inelastic scattering and hadron-hadron collisions have been 
described surprisingly well within this approach. Meanwhile calculations of 
such observables have been carried out to many more quantities (see presen- 
tations in Refs. 15-18. The success of these calculations implies that color 
confinement is governed by rather soft processes which allow the close similar- 
ity of parton and hadron momenta. 

This paper aims at a review of such infrared sensitive quantities for which 
there are analytical calculations. M any of the topics mentioned here are dis- 
cussed in more details elsewhere. H^LLfJ Our main goal here is to survey the 
basic ideas and to illustrate the latest phenomenological advances. So we have 
tried wherever possible to refer to the very recent experimental data. 

The approach discussed here is very restrictive and besides the QCD scale 
A there is only one essential non-perturbative parameter, the effective cut- 
off Qo which is a characteristic of the onset of non-perturbative effects. The 
hope is that the study of infrared sensitive quantities can provide us with the 
important information on the confinement mechanism, besides the test of the 
perturbative QCD dynamics. Furthermore, let us emphasize that a detailed 
understanding of the physics of QCD jets is also important for the design of 
experiments and could provide useful additional tools to study other physics. 
For instance, it could play a valuable role in digging out the signals for new 
physics from the conventional QCD backgrounds. 

2 Perturbative Parton Cascades and Jets 



2.1 Jet Definitions and Multi-Jet Structure 

In high energy hard colhsions the hadronic final state consists of many par- 
ticles. Its characteristic feature is the jet structure with bundles of hadrons 
coUimated along certain (jet) directions. In order to quantitatively describe 
this phenomenon one has to introduce the notion of resolution: the higher the 
resolution of the observation the more jets can be distinguished in the final 
state. 

Resolution dependent exclusive cross sections are already familiar from 
QED. The cross section for all final states which are indistinguishable according 
to a resolution criterion is finite in perturbation theory. This result is expressed 
by the Bloch-Nordsiek theoreniLJ but it has its correspondence in QCD. Finite 
jet cross sections in QCD were studied first by Sterman and Weinberg.ci] They 
counted two partons (jets) as distinguishable objects if their fractional energies 
Xi and relative angles 0^; satisfied the conditions: 

Xk > s, xi > e and Qki > S (1) 

for given resolution parameters e and S. In the case of e~^e~ — > qqg the 
configurations satisfying (|l|) correspond to 3-jet events; since the soft dk/k 
and coUinear dQ/Q singularities from the gluon bremsstrahlung are avoided 
the cross section is finite. In the complementary case the singularities in the 
qqg final states cancel against the singularities from the virtual corrections to 
the qq and the rate for these indistinguishable configurations is finite as well. 
In the last years the resolution criterion based on re]a,tiv e transve rse mo- 
mentum has been widely used {kx or "Durham-algorithm"|) .E^IIrzlErl In this 
scheme two particles (jets) of energies Ek,Ei and relative angle 9^; are con- 
sidered distinguishable if 



yu 



-2{l-coseki)inm{ElEf)/s > y,. (2) 



for a given cut-off parameter y^ where s is the cms energy squared of the 
full event. In this case the coUinear and soft singularities are regularized by a 
single parameter j/c — Q\l ^ which corresponds to the resolution 1/j/c- For small 
angles 0fe; the separation variable corresponds to the transverse momentum of 
the lower energy particle with respect to the higher energy jet particle 

0fei < 1 : Vki ^ k^ki/s > Vc or krM > Qc- (3) 



''This algorithro-was discussed within the framework of analytical calculations at the 
Durham WmpkahopF^. The very idea of k^ clustering has been already applied since the 
early eightieaifa. 



It is an important advantage of this scheme that some observables can be calcu- 
lated in-iall orders of the perturbation series in certain logarithmic approxima- 
tinns p'Jp^ a success not achieved within the earlier JADE algorithm.lr3 Fur- 
thermore, the "hadronization corrections" obtained in popular parton shower 
Monte Carlo's have been found to be quite small in the kx algorithm.E2l 

In the applications considered here we only deal with this fcr algorithm. 
We emphazise the recent development ("Cambridge algorithm" Eil) with a dif- 
ferent treatment of the soft particles and the class of algorithms for ep and 
PP/pP collisions which take into account the spectator jets, discussed in a re- 
cent surveyJ££ 



These resolution criteria form the basis to obtain finite results for jet cross 
sections in perturbation theory. At the same time they can be applied to 
define the jet in the experimental analysis of the multihadron final state by an 
iterative procedure. For every pair of particles one computes the corresponding 
distance yu (restricting here to the case of a single resolution parameter). If 
the smallest distance in the event is smaller than the resolution parameter j/c 
the two particles are combined into a single jet according to a recombination 
scheme. In the simplest prescription (E-scheme) the 4-momenta of jets are 
added; alternatively, one may require that either the momenta or the energies 
are rescaled in such a way that massless jets are obtained. This procedure is 
repeated until all pairs satisfy yui > yc- The remaining objects are the jets at 
resolution l/j/c- 

In this way one can study the multiparticle final state at variable resolu- 
tion between the extreme limits: very narrow jets - ultimately the final state 
hadrons - at high resolution and "fat jets" at low resolution. For example, the 
multiplicity N of jets in e+e^ annihilation approaches 



at high resolution 


{yc- 
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"^ -^*hadrons 


at low resolution 


ivc- 


-1): 


N - 


-> 2 



(4) 

The study of events at variable resolution leads naturally to two different classes 
of observations: 

(i) Multi-jet topologies 

At low resolution (y^ large) there are only a few jets and the cross sections 
can be obtained from low order perturbation theory. In many cases higher 
loop calculations are available. For example, the processes e+e" — > 2 or 3 
jets can be calculated from the relevant matrix elements and they are fully 
known in 0{aj.) which corresponds to inclusion of final states with at most 
four partons; the result will then depend on the coupling constant as, the 



only parameter of the theory, and the resolution parameter j/c which can be 
preset by the experimenter, but should be large enough in order to justify the 
approximation, typically j/c > 0.02. The accuracy of the result and the range 
of applicability in yc can be increased if the corrections from higher orders in 
as are included. 

A key role in particle and jet production in QCD is played by the gluon 
bremsstrahlung. Let us recall first the basic process, the gluon radiation off 
a quark: the differential spectrum of the gluon emitted from the quark of 
momentum p with energy k in an approximation of small transverse momenta 
k± is given by the well known formula 

dw^-^^^ = ^lM2C^ [l + (l-z)2] ^^, 
47r z kj_ 

where Ca = Nc, Cp = [N^ - l)/2Nc = 4/3, iV^ = 3 is the number of colors 
and k^ is the gluon 4-momentum and z = k/p. The strong coupling constant 
Us runs with the scale Q which has to be chosen according to the particular 
problem; A is the QCD-scale and Uf is the number of flavors. In addition, the 
processes g ^ g and g —^ qq contribute to the evolution of the parton cascade. 
A characteristic feature of the bremsstrahlung probability (^) is the broad 
logarithmic distribution over the gluon energy k and transverse momentum k± 
and the soft and coUinear singularities in these variables. Considering multi-jet 
events the relative transverse momenta are large (Qc ~ Q), so the integral over 
the bremsstrahlung spectrum gives a factor of 0(1), furthermore, the coupling 
is small and the perturbation series converges rapidly. 

Topology with extra jet: k± ^ k ^ p -^ prob ^^ — 0{1) <C 1. (6) 

TT 

The rate for multi-jet events will decrease with jet multiplicity like a". 

(ii) Inside jet activity 

If we increase the resolution (lower y^ parameter) more and more jets will 
become resolved with decreasing relative transverse momentum. The integral 
over k± down to a low cut-off Qc will generate large double logarithms from 
the collinear and soft singularities in (ph 

Inside jet activity: k± -^ k -^ p — > prob ~ — In {p/Qo) > 1- (7) 

TT 



Therefore, the probabihty for additional gkion emission is not small because 
of the typically large logarithms. In this case, the higher order terms have 
to be resummed to get reliable results. Besides the q ^ qg process, there 
is another "double logarithmic" splitting g ^ gg and the single-logarithmic 
g —^ qq. These processes build up the parton cascade. 

By lowering the cut-off Qc the number of sub-jets inside a primary "fat" 
jet increases and one may ask how far down in transverse momentum can 
perturbation theory be applied; at some point all hadrons will be resolved and 
the sub-jets coincide with the final hadrons. We will see later on, that indeed 
various features of the hadronic final state can be described by perturbation 
theory for partons with low cut-off Qq of a few hundred MeV. In this review 
our main interest is not so much in the precise description of multi-jet events 
which are important for the accurate determination of as but in the low scale 
phenomena at high resolution inside jets and in between the jets. 

2.2 QCD Cascade and Evolution Equations 

The concept of the evolution equation plays a central role in perturbative 
QCD. Therefore we begin with a short account of how it appeared in the 
context of QCD cascade. After the discovery of Bjorken scalin^rlLill in Deep 
Inelastic Scattering (DIS), and the birth of the parton model,n'Ll the idea of a 
parton carrying a finite fraction x of the parent hadron momentum was well 
established. In the Bjorken limit the structure functions became independent 
of the large momentum transfer Q^ and were directly related to the parton 
distributions q{x), for example 

Quantum Chromodynamics has introduced a subtle modification to the above 
relations. With partons as quarks and gluons possessing nonabelian, asymp- 
totically free interactions, the densities acquired a very gentle, logarithimc 
dependence on the large scale Q involved in a process, q{x) -^ q{x,Q^)E}i 
Subsequent, experimental confirmation of these scaling violations is considered 
as a great thriumph of perturbative QCD. Correspondingly it also marked the 
beginning of further active developments in the theory and its applications. 
Predictions of the scale dependence of the partonic distrib utions are cus- 



tomarily given in the form of the DGLAP evolution equationE^lijtdZ 

^q{x,t)^^ f —P{x/z,as{t))q{z,t), (9) 

at 2tt J^ z 

9 



where t = InQ/fi, with some arbitrary factorization scale /i (see below). The 
splitting function P{x,as) is calculable as a power series in as- Eq. (M is 
generic: the complete theory of gluons coupled to Uf flavors of quarks leads to 
a system of 2nf + l coupled equations with P elevated to the (2n/ + l) x (2n/+l) 
transition matrix. The above equation was derived originally with the aid of 
the operator product expansion and the renormalization group methods in the 
space-like regionJH Subsequently it was established with the diagramm atic tech- 
niques, also in the time-like domain and for higher (two loop) orderj^^L^LlI 
(for a comprehensive account, see for example, Ref. 38. The Q dependence of 
partonic distributions, in a local theory like QCD, is caused by weakly damped 
transverse momenta of partons, c.f. (||). Consequently, they are effectively 
bounded only by the phase space for a given process which introduces the 
logarithmic dependence on the scale Q. Fortunately, all the_kjgadlhnis gen- 
erated in this way obey the powerful factorization theorenLci^'EiilciJ Namely, 
they are independent of the particular hard process under consideration, hence 
they naturally belong to the external partons involved in the reaction. This 
is how the quark densities become Q dependent. The original cross section 
can be split into a universal, scale dependent part, and a non- universal hard 
part which however does not contain large logarithms. The universal part will 
then renormalize the external densities so that the final, renormalized densities 
become dependent on Q. In the process of splitting a supplementary factor- 
ization scale fj, is introduced. The final physical result is independent of /i and 
the Eq. (O) can be considered as a renormalization group equation expressing 
this independence. 

On the other hand, Eq. (|9|) has yet another, physically appealing interpre- 
tation.r^r^ Namely it can be regarded as a master equation for the Markov 
process where, at any moment of a fictitious time t = In Q/ fi, any parton can 
emit another parton with a probability, per unit of time, related to the split- 
ting function P. In this way a QCD cascade develops and the dependence of 
the partonic densities on t is easily understood. Even more importantly, this 
interpretation follows directly from the diagrammatic derivation of Eq. (H). 
In a particular class of gauges, called physical gauges, large logarithms which 
generate the evolution are contained only in diagrams which reduce to absolute 
squares of the amplitudes hence allowing for the above probabilistic interpre- 
tation. Moreover, it was found that even in the case where the interference 
effects are important, giving rise to the angular ordering the net contribution 
can be presented in a probabilistic form (see later). 

The statistical interpretation of Eq. (O) is even more evident upon closer 
inspection of the splitting function P. In the simplest case of the q ^ q 
transition, which is relevant for the evolution of the Non-Singlet densities, P 

10 



reads in the lowest order in a. 



1 



1-x 



P{x)=Cf[- =P(.t)- / Piz)dz\5il~x), (10) 







where P{x) = C_f(1 + x'^)/{l — x) and the last part of the equation is defined 
only in a distribution sense. Hence the Eq. (||) can be rewritten as 



d as(t) 



dz -^ / dy -^ 

—P{x/z)q{z,t)-q{x,t) —P{y/x) 



(11) 



This is the standard form of the gains - losses (or master) equation in stochas- 
tic processes. The change in the number of partons (with given x) is the result 
of two opposite effects: increase due to the emissions from partons with mo- 
menta higher than x, and decrease caused by emissions towards smaller than 



42 



X momentarifl In the original diagrammatic derivation the two terms are given 
by the real and virtual emission diagrams respectively. In fact, each of these 
terms separately is logarithmically divergent corresponding to the infrared di- 
vergence of exclusive (here elastic) cross sections. The final sum in Eq. (^0|) of 
the elastic and inelastic contributions is infrared finiteli^ and is conveniently 
represented by the (. . .)+ prescription. 

Another important ingredient, which emphasizes the probabilistic nature 
of the evolution equation, is the Sudakov form factor S{Q^, fj?) first calculated 
in QEDui and intensively used in perturbative QCD calculations, see for 
instanceli3'ii2l Let us construct a probability that a parton with momentum 
fraction x does not emit any other parton during the evolution betwen ti and 
i2, say. Since P is the probability of emission per unit time, the probability 
that there is no emission in a small interval At reads 

r(t, i + At) = 1 - At^ r ^P{y/x). (12) 

27r Jo ^ 

The probability R{ti,t2) that there is no emission in a finite interval (^1,^2) is 
given by the product 

R{ti,t2)^Y[r{t,,t, + l)^exp(-~ r f ^Mp(^z)dzdt\ (13) 

which is just the Sudakov form factor 

i?(ti,t2)-5(ti,t2). (14) 

11 



Hence the Sudakov form factor, which originaUy was derived by diagrammatic 
calculations, has also the simple probabilistic interpretation. The Sudakov 
form factor is infrared divergent since it describes the elastic process. Therefore 
one usually introduces a small infrared cut off e which defines so called resolved 
emissions. As an example consider a quark with a momentum fraction z > l — e 
emitted from a parent quark. This must be accompanied by an emission of a 
soft gluon with Xg < e. Such a process is considered as unresolved. Hence 

SSiM)=e^v{- j' j '^P{z)dzdt\ (15) 

gives the probability for no resolved emissions in (ti,t2) ^nd consequently 
contains any number of very soft, unresolved gluonsf^ 

Evolution equations can be written in many equivalent forms depending on 
the foreseen application. For example, one can eliminate the virtual emission 



term from Eq. (11) with the aid of the Sudakov form factor 



^ ( +\ "''iW [^ '^^ hi I \ I ^\ , g(^.0 d 



or 



d ( q{x,t) \ _ as{t) [^ dz f,( / ^ <l{z,t) 



f ^p(:,/2)liM 
Jx z S{to,} 



dt\Sito,t)J 2t: J^ z ' ' S{to,t)' 
This can be readily integrated to give 



(17) 



q{x,t)^Sito,t)q{x,to)+ I dt'S{t',t) [ ^^f^P{x/z)q{z,t'). (18) 

Jto Jx Z ^7^ 

This integral form of the evolution equation has again a straightforward prob- 
abilistic interpretation with t' being the instant of the last emission. Note that 
only the real processes, described by P and by the Sudakov form factor, appear 
in this formulation. For that reason Eqs. (^^8|) are the suitable basis of very 
successful Monte Carlo simulations of the QCD cascade. 



iitat 

li 



2.3 Evolution Equations for Jet Observables 



The success of the parton description of Deep Inelastic Scattering prompted 
an avalanche of detailed studies of the QCD cascade. In particular, more and 
more properties of the final states produced in various high energy collisions 

'^In some cases also the very low momentum fractions should be exluded " z > e" if they 
correspond to the unresolved emissions, e.g. in the g ^ gg splitting. 

12 



are being confronted with QCD predictions. This required the development 
of the evolution equations for the time- like region and their generalizations to 
other observables. Complete information about any multiparticle process is 
contained in the generating functionalc3 

Zi{u}) ^J2[ '^^^l • ■ ■ d^knUih) . . . u{kn)Pn{kl, . . . , fc„) (19) 

where P^"'{ki, . . . , fc„) is the probability density for exclusive production of 
particles with 3-momenta ki . . . kn and u{k) is an auxiliary profile function. 
Global characteristics like total multiplicity and its various moments are ob- 
tained by differentiating (n9h with respect to the constant profile parameter 
u. On the other hand the k dependence contains the information about the 
differential densities 

Z?(")(fci,...,A;„) = '5"^(M)/Mfci)-Mfcn) U=i, (20) 

and correlation functions (or cumulants) 

r(") (fci , . . . , A:„) = 5" In Z{{u})/Su{ki) . . . Su{kn) |„=i (21) 

of arbitrary order. 

Originally partonic distributions were calculated in the so called Leading 
Logarithmic Approximation (LLA). It applies for both DIS and e+e^ anni- 
hilaton in the region of finite momentum fractions, 0.1 < x < 1., say. Since 
historically the first calculations were done for the space-like Deep Inelastic 
Scattering we shall explain the basic principle of this approach on this exam- 
ple. Since partonic densities measured there represent the total cross section, 
the infrared divergences cancel. The remaining terms are of the type In x In Q^ . 
In the finite x region they yield single logarithms of the form a" ln((5^//x^)". 
These terms are summed in the Leading Logarithmic Approximmation giving 
for example Eq. (j^ . The In x terms are the remainders of the infrared cancella- 
tions mentioned above. They are small for larger x. Therefore one can also say 
that the LLA is valid when the infrared cancellation between gains and losses 
terms in Eq. (|ll| ) is large, i.e. when the effect of the Sudakov form factor is 
important. As mentioned above the Leading Logarithmic Approximation was 
also extended to the e^e~ annihilation and the evolution equation (pi) can be 
used for fragmentation functions for finite momentum fractions z. The higher 
corrections containing powers of as unbalanced by the large logarithms may 
then be systematically included. 

On the other hand, the more detailed characteristics of the final states, 
like for example the global and differential multiplicities, are dominated by 

13 



the region of small z ^ 1/ ^JQ. This is the region we are mostly concerned 
with in the present review. These soft emissions reveal a new , very character- 
istic feature of the QCD cascade - the angular ordering.EJC^ It is basically 
a nonabelian generalization of the well known ChudakovEil effect in QED - 
the soft radiation from a relativistic e^e~ pair is confined to a cone bounded 
by the electron and positron momenta. Analogously, a soft gluon in a fully 
developed cascade is emitted only inside a cone bounded by the momenta of 
its two immediate predecessors. Mathematically, this is caused by a negative 
interference outside the above cone. Interestingly, these quantum interference 
effects can be again (at least in the large Nc limit) cast into a probabilistic 
scheme described by the evolution equations. 

In the soft region, z ~ l/VQj the Sudakov term in Eq. ( |l6|) is not impor- 
tant in the first approximation. In such a case the infrared logarithms do not 
cancel and each power of as is accompanied by the two (soft and coUinear) 
logarithms, i.e. multiplicities are not infrared safc__T]iea£ contributions are 
summed by the Double Logarithmic Approximation.P lp rrM Together with the 
angular ordering it reproduces qualitatively the most important properties of 
the QCD cascade and retains conceptual and technical simplicity. Hence DLA 
remains an important tool of pcrturbative QCD. 

The Double Logarithmic Approximation is valid quantitatively only at 
asymptotically high, energies. Fortunately an important class of corrections 
which bring the applicability range down to the presently available energies, 
does not spoil the probabilistic interpretation. The new scheme, known as 
the Modified Leading Logarithmic Armroximation (MLLA), contains all next- 



to-leading logarithmic corrections) Ir I and is now considered as a standard 



in quantitative tests of pcrturbative predictions. Formally, it includes consis- 
tently all 0{as) terms in addition to the DLA 0{^/al) contributions in the 
exponential factors. 

The complete MLLA description of the QCD cascade has a form of the 
system of two coupled integral evolution equations for the generating functions 
Za, each describing the cascade oiigiaating from the highly virtual time- like 
parton a = q,g with momentum fl nP^ 

ZaiP,e:{u{k)}) ^ e-^'^P^'^Uaip) 

]_ y- f ^ f ^^g-»a(pe)+^„(pe') 



2! t" Jo e' 







as(fci)r,bc 



Pt{z)Z,{zP, e'; {u}) Z,((l - z)P, e'; {u})e(kl - Ql). (22) 



Simplified versions of this equation had been obtained already before, 
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The Pa'^iz) are the sphtting functionsc^ 

1 






P^il-z)=Tn[z' + il-zn 

1 



P^il-z) = 2CA 



z(l 





(24) 


z 


(25) 



with Cf as in (|) and Tr = 1/2. 

The first term in the r.h.s. of ( p2| ) corresponds to the case when the a- 
jet consists of the parent parton only. The integral term describes the first 
splitting a —f b + c with angle G' between the products. The Sudakov form 
factor guarantees this decay to be the first one: it is the probability to emit 
nothing in the angular interval between O' and O. The two last factors account 
for the further evolution of the produced subjets b and c having smaller energies 
and smaller 6' than the opening angle as required by angular ordering. 

Using the normalization property of the GF 



z,(p,e;M)|„(fc)^i = 1 

the MLLA Sudakov formfactors can be found from 



Wo 



Wo 










Qo/P 



dz 



as{kl) 



Qo/P 



6' 



^.-m 



psi^), 



I p^i-) 



nfPHz) 



(26) 



(27) 



(28) 



Collinear and soft singularities in Eqs. (p2[),(|27|) and (EST) are regularized by 
the transverse momentum restriction 



> Qq 



(29) 
z(l - z)Pe' for 



where Qo is a cut-off parameter in the cascades and fcj 
small angles O'. 

Differentiating the product Za exp [wa{PQ)] with respect to 6 and using 
Eq. ( p^ ) one arrives at the Master Equation 



Itae^"'''-^' 2 



IJlt 



dz 



:^i^P^(z)[Z,(zP,e)Z,((l 



z)P,e) - Za{P,Q)] (30) 
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which gives, as one of the apphcations, the MLLA counterpart of the DIS 
evolution equation (0) for single parton densities. 

Contrary to the differential evolution equations the integral evolution equa- 
tions determine also the initial conditions of the system. In this case they follow 
from Eq. ( p6| ) 

Za{P,Q;{u})\pe=Qo^Ua{P) (31) 

with the simple interpretation that the jet originating from a parton a contains 
only this parton at the lowest virtuality Qo- 

Generating functions Zq and Zg form the building blocks sufficient to de- 
scribe the complete final states realized in physical reactions. For instance, the 
e"^e^ annihilation into hadrons at the total cms energy W = 2P is described 
by 

. ^^--^-^^ - ".qy 

The equations (pG) are now actively exploited and various applications will be 
discussed below. 

In the Double Logarithmic Approximation the parton b, say, emitted at 
the elementary vertex a —^ b + c is considered so soft that it does not influence 
the original parton a. Consequently 1 — z — )■ 1 and c ^ a in Eq.(pfl). This 
soft energy assumption also implies that the splitting functions can be replaced 
by their most singular parts at z ~ 0. This yields a simpler equ ati o n which 
can be integrated to give the DLA Master Evolution EquationElJ'ErfEB (with 
d^k = dLU(fk±) 



Z,+,- {W; M) = [Zq{P, e ^ vr; {u})f . (32) 



Zp(P,e;M) = u{P)exp( J 



duj d'^k± 

LU 27rfc^ 
r(p,e) 

2as{kl) 



[Zg{k,ek;{u})~l]j, (33) 

Cp refers to the respective color factors C'a and Cp- The secondary gluon g 
is emitted into the interval r(P, 9): uj < E — \P\ and 8^ < 6. Due to the 
angular ordering constraint the emission of this gluon is bounded by its angle 
Qk to the primary parton p. As mentioned earlier the DLA, although formally 
correct only at the asymptotically high energies, reproduces satisfactorily the 
general structure of the QCD casacade and allows for an important analytical 
insight. 

At this point we would like to emphasize a deep and beautiful universality 
of the above methods. Even though historically the LLA was first used in 
the Deep Inelastic Scattering and DLA and MLLA in e+e~ annihilation, the 
LLA also applies to the latter at finite z. Similarly the double logarithmic 
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asymptotics also successfully describes the DIS structure functions at large 
In [l/x) ^ In Q^, and MLLA master equation can be used to derive the DGLAP 
evolution equation which is also applicable in DIS. 

At the same time the region lying yet "beyond" the DLA asymptotics in 
DIS, i.e. In (l/x) ;:§> InQ^, is being intensively studied. It is described by 
the BFKL evolution equation and is important for our understanding of the 
emergence Pomeron trajectory in QCD. 

2.4 Parton Hadron Duality Approaches 

At present the application of QCD to multiparticle production is not possible 
without additional assumptions about the hadronization process at large dis- 
tances which is governed by the color-confinement forces. The simplest idea 
is to treat hadronization as long-distance process, involving only small mo- 
mentum tranfers, and to compare directly the perturbative predictions at the 
partonic level with the corresponding measurements at the hadronic level. This 
can be applied at first to the total cross sections, and then to jet production for 
a given resolution; here the partons are compared to hadronic jets at the same 
resolution and kinematics. This approach has led to spectacular successes and 
has built up our present confidence in the correctness of QCD as the theory 
of strong interactions. In such applications the resolution or cut-off scale is 
normally a fixed fraction of the primary energy itself. 

It is then natural to ask whether such a dual correspondence can be carried 
out further to the level of partons and hadrons themselves. The answer is, 
in general, affirmative for "infrared and coUinear safe" observables which do 
not change if a soft particle is added or one particle splits into the coUinear 
particles. Such observables become insensitive to the cut-off Qo for small Qq. 
Quantities of this type are energy fiows and correlations and global event shapes 
like thrust etc. 

In the next step of comparison between partons and hadrons we consider 
observables which count individual particles, for example, particle multiplici- 
ties, inclusive spectra and multiparton correlations. Such observables depend 
explicitly on the cut-off Qq (the smaller the cut-off, the larger the particle 
multiplicity) . 



The very assumption of the hypothesis of Local Parton Hadron Dualitjii 
is that the particle yield is described by a parton cascade where the conversion 
of partons into hadrons occurs at a low virtuality scale, of the order of hadronic 
masses (Qo ~ few hundred MeV), independent of the scale of the primary 
hard process, and involves only low-momentum transfers; it is assumed that 
the results obtained for partons apply to hadrons as well. 
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Within the LPHD approach, PQCD calculations have been carried out in 
the simplest case (at asymptotically high energies) in the Double Logarith- 
mic Approximation or in the Modified Leading Logarithmic Approximation 
which includes higher order terms of relative order yfa^ (e.g. finite energy cor- 
rections); they are essential for quantitative agreement with data at realistic 
energies. According to LPHD, the shape of the so-called "limiting" spectrum 
which is obtained by formally setting Qq = A in the parton evolution equations, 
should be mathematically similar to that of the inclusive hadron distribution. 

In this review we examine, in particular, applications of the LPHD sce- 
nario concerning "infrared sensitive quantities" . To deal with the cut-off Qo 
one can proceed in different ways. If the cut-off dependence factorizes (for 
example, multiplicity) one can again get "infrared safe" predictions after a 
proper normalization. In other cases (for example, inclusive momentum spec- 
tra) the observables become insensitive to the cut-off at very high energies if 
appropriately rescaled quantities are used. 

More generally, one can test parton-hadron duality relations between par- 
tonic and hadronic characteristics of the type of 

0{xi,X2, ■ ■ ■)\hadrons ^ K 0{xi,X2,- ■ ■ , Qo, -^)\partons (34) 

where the non-pcrturbativc cut-off Qq and the "conversion coefficient" K 
should be determined by experiment (for review, see Ref. 17). An essential 
point is that this conversion coefhcient should be a true constant independent 
of the hardness of the underlying process. 

The hypothesis of LPHD lies in the very heart of the analytical perturba- 
tive scenario, but a the same time this key hypothesis could be considered as its 
Achilles heel as it remains outside of what can be derived within the established 
framework of QCD today. One motivation of LPHD is the "pre-confinement" 
property of QCDEZI which ensures that color charges are compensated locally 
and color neutral clusters of limited masses are formed within the perturbative 
cascade. On the other hand, LPHD fits naturally into the space-time picture 
of the hadroformation in QCD jets to be discussed below. 

When comparing differential parton and hadron distributions there can 
be a mismatch near the soft limit caused by the mass effects. This mismatch 
can be avoided bv a proper choice of energy and momentum variables. In a 
simple model23'Ej partons and hadrons are compared at the same energy (or 
transverse mass) using an effective mass Qq for the hadrons, i.e. 



i^T,parton — kT,parton ^ ^TJiadron — \J '^T.hadron '^ ^0 (^5) 

then, the corresponding lower limits are kr.parton -^ Qo and kr^hadron — > 0. 
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Finally, let us recall that within LPHD approach there is no convincing way 
to introduce the different hadron species. For this one must resort to models. 
These are also vital for the practical purposes, for instance, for unfolding parton 
distributions from hadron spectra. At the moment all the so-called WlC ged 
= (With Interfering Gluons) Monte Carlo models (HERWIcil, JETSEtE, 
ARIADNEH-U) are very successful in the representation of the existing data 
and they are intensively used for the predictions of the results of present and 
future measurements. It is worthwhile to mention that for many observables 
the LPHD concept is quantitatively realized within these algorithmic schemes. 

2. 5 Space- Time Picture of Jet Evolution 



To exemplify the space-time structure of the development of the QCD jetsli^nH^Ei! 
let us consider the process e^e^ — > qq. This may be viewed as the decay of 
a highly virtual photon with mass Q, or as a real Z^, the decay time scale 
is short, tannih — l/Q ~ 10~^ — 10~^/m. The qq pair is kicked out of the 
vacuum as bare (at scale 1/Q) color charges until the gluon field has had time 
to regenerate out to a typical hadron size R '^ ifm. Allowing for the Lorentz 
boost this takes a time thad — Q/m x i? « QR^ ~ 10^ fm, where the second 
approximation is appropriate to light hadrons. 

Since thad S> tannih the qucstion arises of how the color charges are con- 
served over the space-like separated distances involved. The primary quarks 
will radiate gluons and here two new scales are relevant. 

First, from the virtuality prior to emission, the formation time of a gluon 
of energy k is tform — k/k"^^. Secondly, for the gluon to reach a transverse 
separation of R and become independent of the emitter takes a time tsep = 
{k±R) ■ tform, whilst the hadronization time may be written as thad{— kR?) = 
{k±)'^ -tform- For such quark-gluon picture to make sense we require k± > R~^ 
so that tform < tsep < thad- Within this scenario the first hadrons arc formed 
at the time t ~ t^rit ~ R- It is the moment when the distance between the 
outgoing q and q approaches R. At t > tcrit the two jets are separated as 
globally blanched, and the parton cascades develop inside each of them. The 
gluon bremsstrahlung becomes intensive only when the transverse distance 
between any two color partons exceeds R. 

With increasing time the partons with larger and larger energies k ^ -^ 
hadronize (inside-outside chain). If fcj_ < R~^ then within perturbative sce- 
nario we can say nothing. On the borderline are quanta with k^R ^ \ (though 
with arbitrary large energies). These do not have enough time to behave as 
free perturbative partons because their hadronization time is comparable with 
the formation time, tform ~ thadr ~ kR^ . 
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We distinguish these quanta from the essentially perturbative gluons by 
the name gluers.LL°l Contrary to conventional QCD partons gluers do not par- 
ticipate in perturbative cascading and their formation is a signal of switching 
on the real strong interactions (a^ ~ 0(1)). 

In this picture soft particles with k ~ R~^ produced at the lower edge of 
the perturbative phase space play a very special role. Their production rate is 
unaffected by the QCD cascading, and in some sense they can be considered 
as the eye-witnesses of the beginning of the "hadronization wave" . 

It is an interesting question, whether the contribution from non-perturbative 
emission (a "pedestal" in the rapidity distribution) not included in the pertur- 
bative calculation (quanta with k± > Qo) manifests itself in the present data. 
The studies_Dn rapidity gaps do not require such an addition to the perturba- 
tive result 1^ within the measurement accuracy aniLsimilar conclusions could 
be drawn from the soft particle production in jets.B^ 

3 Multiplicities 

The simplest global characteristic of the hadronic final state is the particle 
multiplicity. The mean multiplicity of partons in a jet is derived from the 
generating functional according to the general rules rt2G) by single differentia- 
tion. Then one obtains from the master equation (|30|) the following coupled 
system of evolution equations for the multiplicities J\fa in quark and gluon jets 
(a = q,9) 



dNa{Y) 






X [J\fb{Y + In z)+JVc{Y + Ml -z))-J\fa{Y)]. (36) 

The multiplicities J\fa depend on the jet virtuality k w PQ and the cut-off Qo 
or on 

y = l„f A.,„f (37) 

where Q denotes the maximum angle between the outgoing partons b and c 
- for large angles Q the variable Y ^ rj = Iu^k/Qq), k = 2Psin(8/2) have 
also been suggested.!^ The integral boundaries in (|3g) are further reduced 
by the condition k± > Qo- The transverse momentum is usually taken as 
k±_ ~ z(l — z)k, in some applications also k± = min(z, 1 — z)k ("Durham k±^). 
The initial condition for solving this system of equations reads 

J^aiY)\Y=0 = 1 (38) 
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which means there is only one particle in a jet at threshold. 

This set of equations determines the multiplicities of partons in quark and 
gluon jets in absolute terms for a given cut-off parameter Qq. In principle, 
one can solve the equations by iteration starting with (jSg) which yields the 
perturbative expansion of this quantity. In some approximate schemes this 
series can actually be resummed analytically. 

3.1 Asymptotic Behavior 

At high energies the solutions of ( |3q ) can be written 

MgiY)r^cxplj -/{y)dy\ (39) 

where the anomalous dimension 7 has the expansion in 70 ^ y^oj 

7 = 7o(l -ai7o-a27o -«37o ■■•): (40) 

likewise the ratio of gluon and quark multiplicity 

rE.^ = g^(l-ri7o-r27o^-r37o^..). (41) 

The coefficients can be found from (pq) by expanding the A/''s at large Y. The 
leadin^S'lll (DLA) and next to leadingHlll (MLLA) coefficients a^ lead 
to the multiplicity growth expressed in terms of the running coupling 

\njViY)^ci/^^jY)+C2\nas{Y)+C3 (42) 

where C3 is an arbitrary constant and 

ci - x/96^/6, C2 - i + ^nf/b. (43) 

In this next-to-leading high energy approximation (MLLA) Afg ^ Afg since the 
difference would be a y^aj correction as given by the terms ri and 02. The 
asymptotic limit r = Ca/Cf — 9/4 in (ETl) has appeared first in the discussion 
of radiation from color octeL^nd triplet sourcea2_l and in the jet calculus in 
Leading Log ApproximationEB. Terms of higher order have been obtained in 
next-to-leading ordeiEjEil, next-to-next-to-leading ordeilH3'Lll and in one yet 
higher order (3NL0)lLi!l. 

The Eq. (Bq) completely determines the leading term (DLA) and the 



second term (MLLA) in (40) which are important for the high energy behavior. 
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The terms of yet higher order (see also the recent review Ref. 71) are not 
completely determined by the single jet evolution equation ( pq ) because they 
are affected by large angle emission processes. However, these terms are of 
some relevance nevertheless as they include energy conservation in improved 
accuracy. Furthermore, the summation of the full perturbation series allows 
one to take into account the initial conditions ( |38| ) at threshold. These results 
will be discussed next. 

3.2 Full Solution in DLA 

In this approximation only the most singular terms in the splitting functions 
~ 1/z are kept, i.e. P| ~ 2Ca/z, P| =i Cf/CaP^, c.f. Sect 2.3. The recoil 
is neglected, i.e. the incoming parton retains its energy and angle in the final 
state. The maximal angle Q in the parton splitting is then the half open- 
ing angle of the jet. Using the logarithmic variables (pj) and the anomalous 
dimension 

7o(Q)- = ^q7a)' ^--^' b=^Nc--n,, (44) 

we obtain the evolution equation for the gluon multiplicity 

^^^ = £ dy^',{yW,{y), (45) 

where the integral is over the intermediate parton momenta k — zP, or y — 
hi{k<d/Qo)- After second differentiation 

JV;^'iY)-^iiY)N;iY)^0, (46) 

AA,(0) = 1, AA;(0) = 0. (47) 

In case of fixed coupling this leads simply to 

1 fpe^ ^° 



Af,{Y) = cosHjoY) c ^[^j , (48) 

with a power behavior at high energies. This result was found already prior to 
QCD in fixed coupling field theorieaij. For running coupling the solution of 
( p^ ) is found in terms of modified Bessel functions 

MgiY) = l3VYTX{Koif3VX)hipVYTX) + Ioi^VX)Ki{pVX)}. (49) 
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At high energies (large Y) the asymptotic expressions Ii,{z) « e^ /V^ttz and 
Ki,{z) ~ e~^y/Tr/2z apply. Then the second term in ( ^ ) can be neglected 
whereas the first term yields the exponential growth of multiplicity 



AA,(y)~expWl^(y + A), (50) 



corresponding to the first term in ([421). Because the coupling is decreasing with 
increasing energy the multiplicity growth is slower than the power ([48| ) of the 
fixed coupling theory but still larger than a logarithm as in a "flat plateau" 
model. Using this formula for hadrons (Qo ^ C'(A)i A <C 1^) the dependence of 
the cut-off parameter A factorizes and determines the absolute normalization. 
We may also apply Eq. (|4^) to jet multiplicities. This means we consider 
Qo ~^ Qc i^ ~^ ^c) as variable cut-off for the relative transverse momenta 
between jets within the Durham jet algorithm (see Sect. 2.1). At high energies 
with the above approximations we find the behavior in the two limits (Q) for 

at high resolution {Qc -^ Qo) ■ A/" -- (P^Y)^/^ In f ^= ) exp ^/pY (51) 

\f3vKJ 

at low resolution (j/c ^ 1) : JV ^ 1, (52) 

where we used Kq{z) ~ ln(2/z) for small z. At high resolution for Qc -^ A 
(Ac — > 0) the parton multiplicity diverges logarithmically because of the Lan- 
dau pole appearing in the running coupling. The pole is shielded by the cut-off 
Qc = Qo and at this value the particle multiplicity reaches the hadron multi- 
plicity Af/K according to the LPHD prescription ( p^ ) . It should be noted that 
for small A the coupling becomes large ~ 0{1). The higher order terms are 
sufficiently suppressed by phase space so that the perturbative series can be 
resummed as demonstrated for DLA by Eq. (p9|). This overall features of the 
DLA are similar to the more precise calculations and the experimental findings 
to be discussed below. 

3.3 Modified Leading Logarithmic Approximation (MLLA) 

The next order terms are generated by the non-singular terms in the splitting 
functions, the z dependence of the coupling and the inclusion of energy con- 
servation in the parton splitting. One rau rn odify the evolution equation and 
keep only terms of ©(y^ol) and 0(Q;g)l^trn neglecting those of higher order 
assuming as to be small. Then one can derive the multiplicity again from a 

23 



differential equation and express the resultH'lizI in a compact form 

AAg(r, A) = zi (^^ {Ib+i(zi)Kb(z2) + Kb+i{zi)Ib{z2)}, (53) 



16iVc ,,^ , ,, _ /16iVc . „ _ a _ 11 . 2n/ 



This expression preserves the initial condition A/'g(0, A) = 1. A simplification 
occurs in the case A — *■ ( "limiting spectrum" ) where one finds from (B3) using 
Kb{z) -^ V{B){z/2)-^/2 and /^(z) ^ for z ^ the finite limit 

<-(n-r(i?)(|)'"''^'^s+i(^i). (55) 



At high energies this result and in the same way the first term of Eq. (53) 
yield the asymptotic form 



N{Y, A) ^ r-s/2+1/4 exp J^^r (56) 

which is equivalent to (p2). The DLA result is recovered from (m) for B ^ Q. 
On the other hand, within the high energy MLLA approximations assuming 
small as, the logarithmic singularity for A — *■ present in the general equation 
( p6| ) and in (E^) of the DLA has disappeared. Generally, the MLLA results 
are expected to differ from the exact solution in the kinematic region of large 
coupling Ug. A common solution to the MLLA simplified coupled evolution 
equations for quark and gluon jets has been derived as well^iid can be repre- 
sented by expressions in terms of modified Bessel functions.Lj3. 

At threshold the second condition A/''(0, A) = should hold which follows 
directly from the evolution equation (p^) . A shortcoming at first sight of the 
full analytical MLLA solutions is a violation of this threshold condition and 
actually A/'(0, A)' < 0. The reason is again that a^ becomes large near threshold 
and taking only the first two terms of the y/al expansion is not justified. This 
is one of the motivations to solve the evolution equation more precisely using 
numerical methods. 

3.4 Numerical Solutions 

The MLLA approximate solutions have been compared with the numerical 
solution of the coupled system of equations in (|3g).L3 For sma 1 A ~ 0.02 



good agreement is found already shortly above threshold (Fig. ^pli2l): for the 
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Figure Ij— . Comparisoiiliij of different approximations to the master equation for jet 
evolutionO as a function of the jet energy variable Y = ln(K/Q^ (k = 2i?sin(©/2), = 
7r/2); (a) Multiplicity in quark4ets; the exact numerical solutionl'T. the MLLA results from 
the coupled Eqns. for q&cg ietsr I. the solutioiLfor g jets only, multiplied by 4/9 and the linfe 
iting spectrum with A = 0, multiplied, by 4/£(££|; (b) Multiplicity ratio for exact solutionl'^ 



compared isith asymptotie-iralue 9/4221, the MLLA results using-the norm|ali^ation at thresh- 
old (D0E3 and CDFWE3) and asymptotic expansions {GmM and DNE^I) for A = 0.02. 



numerical solution with the exact integration limits the inelastic threshold with 
J\f > 1 starts at y > 0, for the analytic solutions with approximate integration 
limits this occurs already at F = with Af' < 0. A larger discrepancy by a 
factor 2 is found for the "limiting spectrum" with A = 0. 

On the other hand, the ratio of giuon and quark multiplicities is rather 
sensitive to the type of approximation as this difference is a sub-leading effect 
(see Fig. np). At LEP energies {Y ~ 5) the inclusion of higher order terms 
decreases the ratio r from the asymptotic value r = 9/4. The asymptotic 
solutions, "GM"l23and "DN"Ey, which have no normalization condition at any 
finite energy eventually will reach unphysical values r < 1 at low energies. The 
fully resummed results are normalized to r = 1 at threshold. The numerical 
solution takes the lowest value of the ratio r at high energies. 



3.5 Experimental Results on Quark Jets 

Tests of the QCD predictions for multiplicities are available from the final 
states in e"'"e~ annihilations and in the current fragmentation region in deep 
inelastic ep scattering. In MLLA accuracy the asymptotic predictions (^2|) can 
be taken over from the single gluon jet to the e"*"e~ final state. In most fits 
of (B2) to the data the 2-loop formula for as is used although the leading log 
calculations are only accurate to one loop. 
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Figure 2: Mean charge particle multiplicity in e+e annihilation as function of total energy 

compared with the MLLA prediction. The data at energies > Mr-have been corrected for 

increased multiplicity from b and c decays. lj 



The agreement with the data is generahy good. As an example we show 
data in the range from 20 to 180 GeV in Fig. ||lirl The curve represents Eq. 
( p4 ) where N' is multiplied with (1 + dy/a^) to allow for the possibility of 
a next-to-next-to-leading-order correction. Fitted parameters are the overall 
normalization, aJmz) = 0.119 ± 0.003 and d — 1.11 ± 0.39. Recent measure- 
ments at LEP-2lLil follow the extrapolation from lower energies although they 
are systematically a bit at the low side. For example, the fit to data between 
12 and 161 GeV pred icts A/'ch — 27.6 at 189 GeV which is to be compared 
with the measurementlLlI A/'ch = 26.95 ± 0.16(stat.) ± 0.51(syst.). Recently 
a comparison of data with the Drediction including the higher order terms up 
to 3NL0 has been performcdllil Fitting the normalization and the coupling 
leads to a good result with reasonable A in the range 0.12-0.26 GeV. 

Alternatively, one can compare the data with predictions from the full 
parton cascade using the normalization at threshold. The numerical solution 
of the pair of evolution equations (36) with (p8|) has been compared with data 



on e+e annihilation; in this calculation the term of order a^ had been replaced 
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Figure 3: Data on the average jet multiplicity Af at Q = 91 GeV for different resolution 
parameters yc (lower set) and the average hadron multiplicity (assuming Af = fA/'ch) at 
different cms energies between Q = 3 and Q = 91 GeV using Qc = Qo = CL508 GeV in 
the parameter yc (upper set). The curves follow from the evolution equation (pm with A = 
0.5 GeV; the upper curve for hadrons is based on the duality picture (b4) with K = 1 and 
parameter A = ln(Qo/^) = 0.015 (Figure from Ref. 74). 



by the result from the full matrix element of the same order .li^ The predictions 
from this analysis for both, jets and hadrons, are shown in Fig. || together with 
the data. 

The lower set of data refers to the jet multiplicity at variable resolution 
Vcut = [Qc/QY' calculated at fixed total energy Q. For ycut -^ 1 all particles 
are combined into two jets and therefore J\f ^ 2 as discussed in (0). On 
the other hand, for i/cut ~* all hadrons are resolved and J\f -^ Mhad- The 
theoretical prediction describes the data reasonably well down to ycut ^ 10^'' 
and is determined fully by the parameter A. 

An improvement of the description for large ycut is possible if 2-loop results 
are used. For ycut > 0.01 the data are well described by the complete matrix 
element calculations to 0{al). In the region ycut > 10""^ the resummation of 

the higher orders in as is important Mj The calculation at large ycut allows 
the precise determi n a t ion of the coupling or, equivalently, of the QCD scale 
parameter ^-jjg- l^'l^ 

The multiplicity diverges for small cut-off Qcut — > A as in this case the 
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Figure 4: Jet multiplicities extending towards lower ycut parameters; full lines as in Fig. 
for jets, dashed lines the same predictions but shifted ycut — » ycut — Q^ / Q^ ^ccovding to 
the different kinematical boundaries as in ([ial_with parameters as in Fig. H (data from 

opalHeJ). 



coupling as{kT) diverges. The divergence is shielded by the cut-off Qo and 



according to the duality picture in (34) the parton multiplicity represents the 
hadron multiplicity at this scale. It is found from the data that this happens 
for Qc = Qo — 0.5 GeV if the total hadron multiplicity is taken as 3/2 of the 
charged multiplicity. The corresponding calculation at lower cms energies is 
in agreement with the hadron multiplicity data down to Q = 3 GeV with the 
same parameter Qq as shown by the upper set of data and the theoretical curve 
in Fig. 0. Interestingly, the normalization constant in (B4) can be chosen as 



K^l, 



(57) 



whereas in previous approximate calculations using the limiting spectrum (55) 
the value K ^ 2 has been obtained.L3 The parameter K is correlated with 
Qo and can be varied within about 30%. The result K = 1 implies that the 
hadrons, in the duality picture, can be viewed as very narrow jets with low 
resolution parameter Qo ^ a, few 100 MeV.fl 

The behavior of multiplicities in Fig. is close to the qualitative expec- 
tations from the DLA discussed in Sect. 3.2. The running of the coupling 
as is crucial for the results. Namely, for constant Ug both curves for hadrons 
and jets would coincide and follow a power law in the ratio of available scales 
Qcut/Q as in (Eq). With running asikx/A) the absolute scale of Qcut matters: 



'^The precise value of Qo depends on the expression used for the kj_ scale in the argument 
of as and varies typically between 250 and 500 MeV. 
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Q, E GeV 

Figure 5: Average charged particle multiplicity in the current fragmentation region of 

the Breit frame as measured by HI and ZEUS Collaborations at momMrtum transfer Q in 

comparison with e+e^ results at cms energy E* = QE3 



as varies most strongly for Qcut — > A for jets at small ycut and for hadrons 
near the threshold of the process at large ycut (small Q) where again as > 1- 

Recently, data became available at smaller ycut in a wide energy range 
from 35 to 183 GeVl2i!lLiJ and examples are shown in Fig. ^ In the theoretical 
calculation all hadrons are resolved for Qcut —^ Qo whereas in the experimental 
quantities, this happens for Qcut -^ 0. This kinematical mismatch can be 
avoidcdUJ by a shift in ycut according to (^). The shifted (dashed) curves in 
Fig. ^ describe the data rather well whereby the Qq parameter has been taken 
from the fit to the hadron multiplicity before; the predictions fall a bit below 
the data at the lower energies like 35 GeV. The non-pcrturbativc Qo correction 
becomes negligible for Qc > 1.5 

It appears that the final stage of hadronization in the jet evolution can 
be well represented by the parton cascade with small cut-off Qo and with 
the standard 1-loop running coupling. This description clearly goes beyond 
standard perturbation theory which is determined entirely by the QCD scale 
A. The calculation in the soft region rather corresponds to a non-perturbative 
model which involves a hadronization scale Qq. In some kinematic regions 
(small k±) the coupling becomes large as/ir '^ C'(l) but ~ as experienced 
with the analytical DLA results - there is good convergence of the leading log 
summation also in this region as the higher order terms are suppressed by soft 
gluon coherence effects. 
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3. 6 Test of Jet Universality 

Resulta^on quark jet fragmentation are expected to be universal in the parton 
model.a If we consider the particle multiplicity the soft particles are included 
and in this case the reference frame becomes important. In deep inelastic lep- 
ton proton scattering at momentum transfer Q^ the particles in the current 
fragmentation region in the Brcit frame should be compared with the multi- 
plicity of one hemisphere in e+e^ collisions at cms energy E* — Q (see, for 
example, Refs. 1,82-84) Results of such a comparison are shown in Fig. pi 
The DIS results approach those from e+e^ for energies Q > 10 GeV. At lower 
energies processes not available in e^e" annihilation like photon gluon fusion 
are important. The agreement at higher energies confirms the universality of 
the jet fragmentation and the relevance of the Breit frame. 

3. 7 Comparison of Quark and Gluon Jets 

The experimental results on multiplicities in quark jets are derived as half of 
the total multiplicity of e+e^ annihilation or the multiplicity in the current 
region of DIS. These experimental results are well met by the perturbative 
calculations (and their non-perturbative extensions). 

The situation is more dif&cult for gluon jets and the following results have 
been presented: 

1. A fully inclusive measurement is possible at lower energies Q from 
radiative decays of T(IS') which are assumed to proceed in lowest perturbative 
order through T — > ^gg. This yields measurements at Q ~ 5 GeV of the gg 
system.EHI Simi larl y, results have been obtained near 10 GeV from the decay 
T(3S') -^ 7X62^—1 assuming Xb2 -^ gg- At Q ~ 5 GeV the ratio r is still 
compatible with unity. 

2. At higher energies the multiplicities have to be extracted from gluon 
jets in a more complex multi-jet environment, either 3-jet events in e^e~ or 
high pt jets in pp colliders. A simple situation is met again in e+e^ — > 3 jets 
with two quark jets (taken as identified b quark jets) recoiling t oge ther against 
the gluon. This possibility has been pointed out some time agoE3 and worked 
out in detail for e^e^ annihilation.E^ Results on th e m ultiplicity in the gluon 
hemisphere have been obtained in this way by OPALcU and the ratio r is found 
to be r = 1.514 ± 0.019 ± 0.034 at Q - 80 GeV. 

3. Furthermore, results have been obtained from symmetric 3 jet events 
in Y-confimration with a quark in one hemisphere and quark and gluon in the 
other one.^ In this case the gluon jet multiplicity is obtained as difference of 
the 3 jet and the known qq multiplicity as suggested by a perturbative leading 
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order analysis! Ir 'I The results obtained in the intermediate energy region 
interpolate between the Q = 10 and Q = 80 results above. 

4. The ratio r has also been determined from high px jet productioiLwith 
dijet masses in the range 100-600 GeV at the TEVATRON by CDF.II The 
energy dependence of the multiplicity has been assumed to follow the "limiting 
spectrum" formula (pa). Using the known composition of quark and gluon jets 
the energy dependence of the multiplicity is then predicted for a given ratio 
r. The curves in Fig. ||, calculated for fixed ratios, are then compared to the 
data which yields the estimate r — 1.7 ± 0.3 over this energy range. 

The experimental results 1. and 2. do not use further theoretical input and 
have been compared with the asymptotic predictions and with the numerical 
solutions of the evolution equations; both are found not fully satisfactory. 

1. In the first approachlLH the 3NL0 asymptotic expression is fit to the 
gluon jet data at Q ~ 5, 10 and 80 GeV by adjusting the normalization and 
A. Using the corresponding expression for r — Mg/Nq and taking the same 
parameters one can predict the multiplicity in quark jets. This is now found 
lower over the full energy region by 20-30%. Correspondingly the ratio r is 
predicted larger than the measurements. 

2. In the second approachliJ the numerical solution of the evolution equa- 
tions (|6|) is fit to the e+e~ data and the parameters Qo and A are determined. 
The fits are found satisfactory (see Sect. 3.5). The gluon jet multiphcity from 
the same solution are now compared with data. There is a good agreement 
with the OPAL result at Q -- 80 GeV but the prediction for CLEG at Q -- 5 
GeV is too large by 20%. It has been argued that at low energy the exact 
treatment of the 0{as) term is important as was found explicitly for quark 
jets. For the gluon jets at the T such calculations have not been performed 
yet. Another explanation could be a large non-perturbative effect for gluons 
which "freezes" the gluon degrees of freedom at low energies. It is interest- 
ing to note that the result from the HERWIG Monte Carlo shows a similar 
behavior with too large ratio r at low energies. 

At higher energies the numerical solutions yield a rise from r ~ 1.54 at 
Q = 100 GeV to r = 1.68 at Q = 600 GeV which is consistent with the CDF 
result within the large error. 

Taking these results together there is a clear evidence for the difference of 
quark and gluon jet multiplicities. The ratio r at presently accessible energies is 
much smaller than the asymptotic value r = 9/4. It is a considerable success 
that at the higher LEP energies there is agreement with perturbative QCD 
calculations if all non-leading logarithmic terms from the MLLA evolution 
equation are included. 
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MLLA predictions for various r 

O CDF Data, cone 0.466 
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KLpHiy= 0.55±0.06±0.08 
r= 1.7±0.3 
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Figure 6: Charged multiplicity in a high pT jet with cone opening 0.466 rad for different dijet 

masses in comparison with MLLA predictions for different values of the ijatio r = Afg jMq ■ 

The curves corresponds to r = 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.25.E3 



4 Momentum Spectra of Particles Inside Jets 

We start from the MLLA evolution equation for particle spectra and consider 
solutions relevant for the small x region. The early DLA results on the ap- 
proximately Gaussian shape of the distribution ( "hump-backed plateau" )EilfBiJ 

in the variable 

1 p 



e = ln. 



E. 



(58) 



'ie* 



for momentum p of the particle has been an important milestone in the applica- 
tion of perturbative QCD to the multiparticle phenomena. They demonstrated 
the relevance of angular ordering and color coherence to the particle spectra 
at low momentum. The success of the improved calculations within MLLA 
accuracy gave strong support to the concept of LPHD. This Gaussian shape 
is found quite universally in jets from various collision processes and also for 
all particle species; however, the differences observed for various species is not 
yet quantitatively understood. We will discuss here some MLLA based com- 
putations and compare them with typical recent results from e"'"e~, ep and 
pp colliders whereby we include some striking phenomena with the very soft 
particles (< 1 GeV), as well as the energy evolution of spectra from threshold 
to asymptotic energies. 

Whereas this analysis is focused on the small x domain of the particle 
spectra there is a complementary approach which treats the evolution equa- 
tion in the larger x-p region (say, x > 0.1) ("fragmentation function") for which 
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results in 2-loop accuracy for the running coupling are available. Such sti 
are an important tool for the determination of the strong couplingH3'c 
This approach is expected to be more precise at large Xp as compared to the 
1-loop MLLA calculations, on the other hand it does not take into full account 
the contributions important for the soft region. Furthermore, in this approach 
an ansatz involving several parameters is required for the fragmentation func- 
tion(s) at a particular finite energy whereas the MLLA/LPHD approach has 
only one essential non-perturbative parameter Qq. 

4-1 Evolution Equation for ^-Spectra and Approximate Solutions 

The inclusive distribution of partons 6 in a jet from parton a is obtained from 
the generating functional through 

xDlix) = E,-^Za{Ea,Q; {4)U=i (59) 

Su{kb) 

and its evolution equation from Eq.(pQ). In the approximation 1 — z « 1 it 
takes the form 

±xD^A-,Y)= Y: f'^'^^'^Aiz)[-^DE{-^,Y + lnz)], (60) 

C=q,q,g ° 

where Y = ln(F8/(3o) and $2 stands for the DGLAP splitting functions. The 
boundary condition for (pOh reads 

xD^{x)\Y=o = S{l-x)5'i. (61) 

The scale of the coupling is given by the transverse momentum taken as k^_ ~ 
z(l — z)EQ. The shower evolution is cut off by the parameter Qo, such that 
fcj_ > Qo and this restriction is understood in (pO). 

The integral equation can be solved by Mellin transform 



D^{Y, A) = / —x^ [xD{x, y. A)] = / d^e-^^ D{£,, Y, A). 
Jo 2; Jo 



(62) 



In flavor space the valence quark and (±) mixtures of sea quarks and 
gluons evolve independently with different "eigenfrequencies" i'±{ijj). At high 
energies and x <^ 1, the dominant contribution to the inclusive spectrum 
comes from the "plus" -term, which we denote by D^{Y, A) = D^{Y^ A). In an 
approximation where only the leading singularity plus a constant term is kept 
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i^+{uj) ~ 4:Nr: /lO" a, one obtains in a sequence of steps from ( |60| ) the following 
evolution equationlL2l 

+^)^^-(^'^)-^^-£^-(^'^)--'^(-+^)£^-(^'^)' 

(63) 
where 

<. = f«c + ||. (04) 

and the boundary condition corresponding to ( pi] ) is D^{0, A) = 1. Restricting 
to the leading term by setting a — 0, i.e. dropping the r.h.s. in (p3), yields 
the DLA evolution equation. One can now derive asymptotic solutions at high 
energies in leading-order and next-to-leading order or a full solution of the 
evolution equation (p3h including the boundary condition at threshold. 
Introducing the anomalous dimension 7^^ according to 




D^iY, A) = D^iYo, A) exp / dy^^ [a^iy]] , (65) 



the evolution equation ( |63| ) can be expressed in terms of a differential equation 
for j^ 

[uj + 7^)7^ = -fi(as)- — 7^ - a(w -t- 7c^) — , (66) 

Ztt das ^TT 

where l3{as) = -^ as(Y) ~ — 6 ^. For the DLA one finds 



(a; + 7^7^ - 7o' = 0, 7.:/""(«s) = 2 (^-^ + V^' + 4^0 J ; (6^) 

with 

-fl^-fl{as) = 2Nc — . (68) 

TT 

Choosing the (-I-) sign for the square root yields the solution which dominates 
at high energies whereas the solution with the (-) sign would die out. 

The next-to-leading order result in MLLA follows from Eq. (^6|) including 
the r.h.s. where the first term proportional to the /^-function keeps trace of 
the running coupling effects while the second accounts for the hard corrections 
to the soft singularities in the splitting function. In comparison to the leading 
term they are of relative order ^/al and the MLLA correction to DLA reads 



DLA , f^ 



ll 
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0[al'^). (69) 



Instead of these asymptotic solutions one can also find directly an analyt- 
ical expression for the exact solution of the differential equation (E^) in terms 
of confluent hypergeometric functions LLj'Ell 

DUY,X) = ^^^^zizim-A + B + l,B + 2,-zi)^{A,B+l,Z2) 

+ e^^-^^ {B + 1)^{A + l,B + 2, zi)<i>{-A + B + l,B + l, -za)}, 

(70) 

where we have used the notation 

ANr a 

A=— ^, B = -, zi=Lu{Y + X), Z2 = iuX. (71) 

buj 

From these representations of D^ one can reconstruct the x (or ^) distributions 
by inverse Mellin transformation 



[xDix,Y,X)]=Di^,Y,X)= f 

J e — i 



— x-i?(c.,y,A) (72) 

-too 27ri 

where the integral runs parallel to the imaginary axis to the right of all singu- 
larities of the integrand in the complex cj-plane. 

4-2 Moments 

It is convenient to analyze the pro perties of the ^ spectrum in greater detail 
in terms of the normalized momentsMBJril 



1 



?, ^< C>=T-r I d^^'DiO- (73) 



Also one defines the cumulant moments Kq or the reduced cumulants kq = 
Kq/a'^ , which are given for 9 < 4 by 

K, = 1= i^ 

K2 = o" =<(e-q'>, 

Ki = ka^ =<{^~-lf> --ia^, (74) 

where the third and forth reduced cumulant moments are the skewness s and 
the kurtosis k of the distribution. 
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The cuniulant moments can be found using the expansion of the Melhn- 
transformed spectrum Di^{Y,X) in (p2 



lnDUY,X) = J2 ^iiY,X) 






(75) 



The high energy behavior of moments in next-to-leading order can be obtained 
with Kg{Y, A) = {—d/dioY ^^D^{Y, A)|„ „■ From the high energy approxima- 
tion of the anomalous dimension 7^^ in (|69| ) using (|65| ) 



K„ 



dy [-Q^] I'^i'^siv)) 



(76) 



L'^O 



which shows the direct dependence of the moments on as{Y). For fixed ofj, for 
example, one obtains simply Kq{Y) ex Y for high energies. Alternatively, one 
can derive evolution equation for the momentsMSI and derive the high energy 
behavior in next-to-leading order. Approximate forms for the yet higher order 
terms can be obtained from Eq. ([70|) for arbitrary A and in particular for the 
limiting spectrum with A = 0. 

The mean of thp ^Histrihutinn in MLLA proves to have an energy depen- 
dence of the forml^"lr ' tr ^ 



^ = Y 



with 



C 



1 C a - .,,„ 

2 V y y 



= 0.2915(0.3513) for 



3(5). 



(77) 



(78) 



16iVc6 

The moments evolve in different energy regions according to the relevant num- 
ber of flavors n/. It turns out that even at the energies of LEP-2 the approxi 
mation nj = 3 is rather good, in any case better than Uf = 5.L3I2M 



A quantity closely related to ^ is the position of the maximum ^* which 
differs fraiii^£_QiiIy in higher order terms a, h in ( [77| ) . For the limiting spectrum 
one findsl 



r 



Y 



C 
Y 



(79) 



This form leads to a nearly linear dependence of ^* on y. It is worthwhile to 
mention that in the large Nc limit, when lliVc ^ 2n/ (cf. Eqs.(||) and (p4|)) 
the parameter C becomes independent on both n/ and Nc and approaches 
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its asymptotic value oi C — ^^ ~ 0.23. Therefore in this hmit the effective 
gradient of the straight hne is determined by such a fundamental parameter of 
QCD as the celebrated ^ factor (characterizing the gluon self interaction) in 
the coefficient b. 

The shape parameters a,k, s for the limiting spectrum have been derived 
in next-to-leading order at Nc = 3 as 

^ ^l- f] + 0{Y-'/^) (80) 



Viz \ 4:Z 

16 a 

where z — ^IdNcY/b. It is worthwhile to notice that the next-to-leading 
effects are very substantial at present energies. In particular, the spectrum 
significantly softens because of energy conservation effects. This influences 
the rate of particle multiplication which is strongly overestimated by the DL 
approximation. 

The asymptotic behavior (n > 1) can be obtained from ( [Tq ) with ( |69| ) 

a^ - y^/2 (83) 



k2n+2 - [VY] (84) 



^2n+l 



n-1/2 

Y] . (85) 



One concludes that the higher cumulants {n > 4) appear to be less significant 
for the shape of the spectrum in the hump region S ^ 1. 

The higher order terms in the series expansion (KQ-R2|) left out are still 
numerically sizable at LEP-1 energies U3 {'^ 10% contribution from next-to- 
MLLA corrections to ^ and cr^) and increase towards lower energies. Therefore, 
it is appropriate in a comparison with the data over a larger energy interval 
to use the full result from the MLLA solution ( |7(]| ) including the boundary 
condition (^l]). A further discussion of the higher moments is found in the 
review .llll 

4-3 Predictions for the ^-Spectra 

Next we present analytical results for some interesting limits, 
(i) Asymptotic Gaussian 
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The simplest example is the DLA prediction which corresponds to the spectrum 

at verv 

shapeE§l^ 



1 energies. Near the maximum one finds an approximately Gaussian 



where 



S 



(e-O 



e 



V2' 
Y 



: exp 



TTCT^ 



-¥' 



{¥ + A)^/^ - A3/2 



(86) 



(87) 



with multiphcity A/'(F) ~ exp(^16A'^c'(^ + A)/6) at high energies as in (pO|). 
This approximately Gaussian shape ( "hump-backed plateau" ) is an important 
prediction of QCD. The drop of the spectrum towards small momenta (large 
^) is a consequence of the coherent emission of soft gluons from the faster ones 
in the jet and we will come back to this phenomenon in more detail below. 

(ii) Limiting Spectrum 

In perturbation theory one is usually restricted to regions where as is small, 
this would require A ^ 1. However, one can see that in Eq. (|8q ) the shape has 
a smooth limit for A — > 0, i.e. as — *■ c». Therefore, the shape is in this sense 
infrared safe. 

In this limit A — > one can derive an analytic expression for the ^ spectrum 
from the full MLLA equation (|7(]|1 using an integral representation for the 



hypergeometric function $,1 



D'"^ii,Y) = 



4iVc 



Xlf_ 



' TT 



-Bo 



cosh a + (1 — 2Q sin h a 



-1 B/2 



4Nc Y g 

b sinh 



'16iVc 



Y 



sinh 



[cosh a + (1 — 2^) sinh a] 



Here a = ao + i£ and ao is determined by tanh ao = 2^ — 1 with C, = 1 — -p. 
Ib is the modified Bessel function of order B. In the present approximation 
this distribution describes the gluon spectrum in a gluon jet. For the quark 
jet this distribution is to be multiplied by Cp/Nc — 4/9. The spectrum (88) 
reproduces the Gaussian behavior near the maximum. 

There is one caveat on the distribution (p8|). As pointed out in Sect. 3.2 
the multiplicity M ^ J D{£^)d^ diverges logarithmically for A ^ 0. In the 
MLLA evolution equation the terms beyond next-to-leading order in the y/a^ 
expansions are dropped which is justified for small as- The gluon emission 
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Figure 7: Spectrum in § = ln(l/a;) as obtained from the MC calculation for a parton 

cascade including the soft gluon coherence which shows the approximatelw-Gaussian shape. 

Also shown is the spectrum for a cascade without coherence .liiliJ 



near threshold a,t ^ ^ Y (small monientuni p ^ Qo) involves large as and is 
therefore not expected to be well approximated in 

(iii) Distorted Gaussian 

The spectrum near the maximum can be represented by a distorted GaussL 
distribution and one finds for small (5 < 1 in terms of the cumulant momentsi 



D{0 



crV27r 



exp 



V ' Q 24 



'-s5 



kS* 



(89) 



Relative to the leading-order predictions, the moments in MLLA accuracy (pQ)- 
(p3) imply that the peak in the ^-distribution is shifted up (i.e. to lower a;), 
narrowed, skewed towards lower ^, and flattened, with tails that fall off more 
rapidly than a Gaussian. 

(iv) Numerical results 
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The evolution of the parton jet in its probabihstic approximation but taking 
into account soft giuon interference can also be treated numerically with Monte 
Carlo methods. In Fig. we show the ^ spectrum for e~^e~ annihilation at 



energy Q — 100 GeV as obtained from the Marchesini-WebberliiL!J MC which 
includes coherence in the soft gluon radiation. One can clearly see the typically 
Gaussian shape. On the other hand, if the coherence is "switched off" in the 
MC the distribution looses its Gaussian shape and instead peaks near the edge 
of phase space for the soft particles {^ <Y, p > Qq). 

(v) The Soft Limit of the Particle Spectrum 

In this limit the coherence of the soft gluon emission plays an important role. 
If a soft gluon is emitted from a qq two-jet system then it cannot resolve with 
its large wave length all individual partons but only "sees" the total charge 
of the primary partons qq. In the analytical treatment, this property follows 
from the dominance of the Born term of 0[as) and one expects a nearly energy 
independent spectrumlij 

This property has been studied recently in greater detail.EaEM An analyt- 
ical formula applicable for the low momenta p and px has been given and its 
general behaviour reads 

for rapidity y and momentum fraction x where the second term is known within 
MLLA and vanishes for pT — > Qo- In this way the energy independence of the 
spectrum follows in the soft limit. 



A crucial predictionErl from this approach is the dependence of the soft 



particle density on the color of the primary partons in (Bffl) : In the soft limit 
the particle density in gluon and quark jets should approach the ratio 

i?(5/9) - §J - ^ for PT^Qo (91) 



i.e. for the minimal px or p. As discussed above, the prediction (91) has been 
obtained in the DLA for the total multiplicities at asymptotic energies. For the 
soft particles it is expected already at finite energies! The consequences have 
been worked out also for soft particle production in DIS and hadron hadron 



collisions. ErI Furt 



recent review 
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ler results on the soft particle production are discussed in a 
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Figure 8: Distribution in fp = ln{l/xp) for charged particles in e+e — > hadrons. The 

full lines represent the result of the simultaneous distorted Gaussian fit to all encrgios-uaing 

moments within MLLA accuracjCj as performed by the DELPHI Collaboration J^^'^ 



4-4 Comparison with Experimental Data 

We begin with a discussion of the ^ spectra. The observation of the Gaussian 
shape of these spectra and a good agreement with MLLA predictions in the 
PETRA energy range was a first success of the LPHD concept .lia The approx- 
imately Gaussian shape is confirmed in the meantime for the spectra in jets 
in the full variety of the hard collisions studied and also for different particle 
species. In general, the limiting spectrum provides a fairly good description of 
the data. We emphasize a few recent results: 

In e^e~ annihilations the study of spectra has been cont iuued towards 
the higher energies of LEP-2. The DELPHI CoUaborationli*^ has fitted the 
distorted Gaussi an paran retrization (p9) to their data up to 183 GeV and those 
at lower energicaiil3'li_2l down to 14 GeV using the moments in next-to-leading 
order ( [77| , |80| -p2|) . The fit has been performed to the central region where the 
distribution was larger than 60% of maximum height. The fitted parameters 
are A — 210 ± 8 MeV, the constant a in ( |77| ) and the normalization using 
71/ = 3. As can be seen from Fig. || the shape of the hump at the different 



41 






Mjj=378 GeV 

6cone=0.28, 0.36 and 0.47 rad 



k 










^^=log4) 



Figure 9: Inclusive g spectrum for charged particles in different jet opening cones fitted 
to the limiting spectrum in MLLA at the dijet maas-of 378 GeV as measured by the CDF 

collaboration.rn 



energies is rather well described. 

The limiting spectrum typically fits i n a b roader ^ range an d gives a good 

the recent re- 
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alsotoLEP-1. 



description of the spectra up to the LEP-lLI 

suit from OPALLiI at 189 GeV indicates a slightly broader distribution (about 

10% larger width) than expected from the extrapolation of lower energies. 

Jets with high transverse energies Et from 40 to 300 GeV in pp colli- 
sions have been investigated by the CDF CoUaborationEHl at the TEVATRON 
whereby also the jet opening angle 8 has been varied between 8 = 0.168 and 
8 = 0.466. The fitted value for A = Qo hi the limiting spectrum was found 
rather stable Qq = 240 ± 40 MeV for the 45 data sets (5 cones x 9 dijet mass 
cuts). As an example, the evolution of the spectrum with the jet opening angle 
is shown in Fig. 0. The curves are seen to represent the data well within a few 
%. Similar results are obtained at other energies whereby the ^ range of the 
successful fit increases with jet energy and opening angle. 

A comparison of the peak position ^* for the various jet energies and 
opening angles is shown in Fig. HOl The spectra depend on the evolution 
variable Y = ln(P8/Qo) for small angle; a comparison is also performed with 
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data at full angle ~ 7r/2 from e+e^ and ep collisions whereby the variable 
Y — ln(Psin Q/Qq) has been used. The data scatter around the expected curve 
( [79| ) for Uf = 3. Taking instead the scaling variable Y = ln(2Psin(8/2)/Qo) 
the full angle data would be shifted to the right by a factor 2sin(7r/4) ^ 1.4. 
This would correspond essentially to a change of the next-to-next-to-leading 
order term in ^79) but would not change the slope. 

The slope is nicely confirmed and the leading DLA contribution (^* ^ Y/2) 
is shown for comparison as the lower curve in Fig. |lO|, adjusted in height; the 
upper curve represents the spectrum for the incoherent cascade which peaks 
near the maximum (^* ^ Y). Apparently the data support the prediction 
from the parton cascade with suppression of soft particles due to coherent 
gluon emission in a large energy range 2Ejet sin ~ 4 — 300 GeV. 



The analytical results for the particle spectrum near the soft limit (90) 
are nicely confirmed by the data. In these calculations the model ( |35| ) for 
mass effects has been used. The experimental data from the available range 
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of energies in e+e annihi lation (Q = 3 — 183 GeV) and Deep Inelastic Scat- 



tering (Q < 20 GeV)E3Ejll23't23 are well described by the calculations. The 
data are consistent with approaching the energy independent limit for small 
momenta p —^ and this conclusion does not depend on choosing a particular 
model for the mass effects. 



It is more difficult to test prediction (|9l|) as a gg system is not easily 
available. From the study of 3-jet events in e+e^ annihilation with the gluon 
recoiling against a qq jet pair an estimate of th is ratio of the soft particle yield 



in gluon and quark jets can be obtained.lLil"! It approaches R{g/q) '^ 1.8 for 
p < 1 GeV which is above the global multiplicity ratio '^ 1.5 in the quark and 
gluon jets but still below the ratio Ca/Cf = 9/4. It is plausible that this 
difference is due to th e fij iite angle between the quark jets, i.e. a deviation 
from exact collinearityJiild It will be interesting to study the radiation pattern 
in more detail. A first result from such studies will be discussed in Section 7. 



5 Multiparton Correlations Inside Jets 

5.1 The General Structure of n-Parton Angular Correlations 

The angular ordering is one of the most characteristic features of the parton 
cascade, c.f. Sub-section 2.3. This basic property of the elementary emissions 
has a direct counterpart in the angular dependence of multiparton correlations. 
They were derived in the analytical form by solving corresponding evolution 
equations with the aid of the resolvent (Greens function) method. This led 
to the integral representations used also in the jet calculus approach.l22l The 
remaining integrals can be solved in the limit of the strong angular ordering 
providing a closed form of many multiparton distributions. For example, the 
fully differen tial angular correlation function of two gluons in a jet (P, 8) 
reads 



iogJirTol 



r(2)(rji,r!2) = D''^\ni,n2) ~ D^^'Hni)D{n2 
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1 f^l 
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^ il^j e(*. -.„)), (92) 

where i^i denotes the angle between the momentum of i-th final parton and the 
momentum of the initial parton P and E = \P\. The two terms in ( P2| ) describe 
two possible ways to produce the final configuration, namely P ^ 1 ^ 2 and 
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P ^r 2 ^ 1. The relative angle 'i?i2 is always restricted by the previous emission 
angle, i.e. dij < "di^i = 1,2, differently in both terms. The above result was 
derived as the high energy limit of the DLA with the constant a^. However 



the generic structure of Eq. (92) remains the same also in more quantitative 
approximations and for arbitrary number of partons. This is seen in the general 
DLA result for n-parton connected correlation function also in the running Us 
case 



r(")(M)c.(^i-)"if^2ii|!!!lexp(^2/?^(e.,n)VlnSVA)Ff({x}). 

(93) 
Again, all the relative angles are limited by the corresponding polar angles 
of previous emissions. Here e^ — \n['di/{-dij)min)/^'^{E'di/ h) and F^{{x}) is 
homogeneous of degree p = —2{n — 1) in the relative angles built from factors 
1/^f,. Similarly to the two parton case, Eq. (|9^), the correlation is a sum of n 
terms, the i-th term describing a configuration with the parton i emitted from 
(connected to) the initial parton at the polar angle "di. All other partons in 
such a configuration are either connected to i (factor i?~- ) or among themselves 
(factor "(9"^ ). The function w(e, n) is a particular feature of the running a^ and 
will be discussed later. For small e (large angles) the simple power behavior 
analogous to Eq. (p3) emerges. 

The above power dependence of the connected correlations is another char- 
acteristic property of a QCD cascade. From a simple picture of a branching 
process one may expect that the resulting parton distributions have a self sim- 
ilar, fractal structure. Indeed, Eq. (02) proves that such expectations are 
correct — QCD cascade is strictly self similar for constant as case. The evo- 
lution of a QCD cascade can be also easily visualized in a fractal p hase space 



119 
which has been constructed within the LUND dipole cascade modelliifl. The 

Renyi dimension\j of the inclusive n-parton configuration has been derived in 

Refs. 109,111,113,114 as 

n -\- 1 
Dn = ^70. (94) 

n 

On the other hand, the running as introduces a scale in the problem, hence 
the self similarity is only approximate. 

Dynamics of QCD provides then a unique prediction for the structure of 
the multiparton correlation. This distinguishes between various phenomeno- 



^This generalization of the HausdorfT dimension allows the fractal dimension to depend 
on n. The Hausdorff dimension reduces to the integer, geometrical dimension for non-fractal 
objects. 
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logical modelsli_L2l of multiparticle production. In fact the way how the n-body 



correlation emerges from the |,wn- hody ones is very similar to the model of Van 
Hove proposed ten years ago 
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5.2 Distribution over the Relative Angle - Comparison with the Experiment 



Predictions of the type of Eq. (|9^) allow to test the Local Parton Hadron Du- 
ality on a more differential, hence more detailed level. However full multibody 
correlations of higher order are notoriously difhcult to measure experimentally^ 
Nevertheless the measurements of particular sections of the two-body corre- 
lations inside a jet already exist and provide important checks of the theory. 
The distribution over the relative angle "(912 is such an object, hence we will 
discuss it in more detail. Consider two partons inside a jet {P,Q),E = \P\ 
at fixed relative angle "(912. The density of z9i2 has been calculated in various 
approximations. In the simplest case (same as in Eq.(|9^)) it reads 



d-di2 

2^12 V Qo y Ul2, 



(95) 



The angular ordering condition takes here the form Qa/E ^ ■!9i2 < 8. A 
more complex expression valid under a weaker restriction Qo/E < "(912 < 6 
is also available. Again the fractal nature of the well developed QCD cascade 
manifests itself via a simple power dependence on the relative angle. The 



absolute density, Eq. (95), depends on the cut-off Qo- Therefore it is advisable 



to define normalized, Qo independent distributions 

diV(^)/d^l2 .., . dM^yd^^2 ,.., 

r[^i2) = T^T— , ri-du) = -^ 96) 

In the first ratio, r, the normalizing quantity is defined as in ( |95|) but with 
D^^^{ni,n2) replaced by D^^H^i) £''^H^2)f]in the second observable, f, the 
square of the particle multiplicity N'^{E, Q) in the forward cone is used as the 
normalizing factor. DLA predictions with the running a^ for these quantities 
read 

r(i9i2) = exp(6(cc;(e,2)-2VT^7)), (97) 



^On the other hand, various higher moments of multiphcity distributions have been mea- 
sured. We do not review this subject here. 

** In the experiment one takes the angle i?i2 between particles of different events ("event 
mixing" ) . 
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r(e) = feexp(6(tj(e,2)-2)), b = 2f3y^\n{EQ/A) (98) 

with /3 as in Eq. (E4). Therefore the dependence on the relative angle enters 
only through a new scaling variable 

ln(eM2) h.{Ee/Qo) 

In {EQ/A) ' - - In (SG/A) ' ^ ^ 

which uniquely combines the three experimentally controlled parameters E, Q 
and "(512. The function uj( e. n] is a solution of a simple algebraic equation 



hence its behavior is known.tiiJ (to good approximation uj{e, n) = n^/1 — e(l - 

(2n^)^^ In (1 — e)). In particular, for the developed cascade, i.e. for '1)12 ~ &, 

e ~ 0, w is linear 

n^ — 1 

uj{e,n)^n~— e, (100) 

2n 

and we recover the power scaling for normalized correlations 

which coincides with the constant a^ result provided one replaces the fixed 
Us by the running as at the scale Q = EQ/ A. 

Both ratios, Eq. (Dq), have been measured no w by DEL PHI in e+e~ and 
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by ZEUS in the current fragmentation region of DIS.ll_LlJ't_L£l The first quantity, 

^(^12), is more sensitive to the genuine correlations as it measures the devia- 
tions from the distribution of uncorrelated pairs, but it depends more critically 
on the choice of the jet axis in the definition of the angles 51^. The second quan- 
tity, r(i?i2), depends only weakly on the jet axis through the opening angle 0, 
however a large (but not a whole) part of its 1^12 dependence follows from the 
simple kinematics of uncorrelated pairs. Experimentally, finding a jet axis is 
not a problem for ZEUS where it is determined by the virtual photon. In 
DELPHI analysis the sphericity axis was used. Scaling in e, revealed in Eqs. 
( P7| , P8|) , is a strongly constraining feature of QCD predictions. In particu- 
lar, it can be tested by changing both Q and E at fixed e. This is suitable 
for DIS kinematics where the momentum of a recoiled jet in the Breit frame, 
P = Q/2, can be readily controlled. Hence the ZEUS data span a range of 
jet energies £' < 30 GeV and DELPHI provides results for the two highest 
energies _E = 45.5 and i? = 91.5 GeV. Both groups use 15° < 6 < 90°. 

ZEUS results for r(i?i2) in a scaling form are displayed in the left columns 
of Fig. n^. Indeed the data show only a residual dependence on E at fixed e. 
The shape of the scaling curve is reproduced rather satisfactorily for higher jet 
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energies. Data clearly indicate the running of as- It is worth noting that the 
running observed here is in the energy range ~ 1 GeV which is substantially 
lower than in other applications. Yet the perturbative description is satisfac- 
tory. DELPHI results (cf. Fig. [l2|a) confirm the QCD scaling predictions 
on f('!?i2) at their highest energies as well. Both groups report that the 8 
dependence is not strong, nevertheless the e is favored as a scaling variable. 

The situation is more complex with the genuine correlations r(t?i2). ZEUS 
and DELPHI results are shown in the right columns of Fig. O and Fig. 12 



respectively. Clearly at low energies ZEUS data have no resemblance to the 
high energy predictions. However a clear transition to the scaling form is seen 
at higher energies. 

DELPHI high energy data ( between E = 45.5 and 91.5 GeV ) indicate 
a scaling behavior for large angles (e < 0.25) albeit below the DLA predic- 
tion. Considerable scaling violations persist throughout the full energy range of 
ZEUS and DELPHI at small angles (large e) . Apparently, the scale breaking is 
related to the normalizing N^J^^^ in ( p6| ) . This factor is built from single particle 
densities D{-d) and strongly depends on the jet axis. This dependenc e can be 
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reduced by measuring the energy-multiplicity-multiplicity correlations 
(see next section). 

Also the theoretical part is based on the rather crude Double Logarithmic 
Approximation. Complete phase space, full matrix element and next to leading 
corrections are not included in (p9). 

With this in mind we find the present state of the confrontation between 
theory and experiment rather satisfactory. At the same time we expect that 
it will stimulate further effort of both parties. 
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5.3 Azim,uthal Correlations to the Next-to- Leading Order 

Complete next-to-leading results are available only for the azimuthal 
tions among gluons emitted from a hard qq pair in an e+e~ annihilation 
This problem is closely related to the string-drag effect and will be discussed 
in detail later. In short, consider emission of a soft gluon associated with the 
hard qqg system. Particle density is reduced around the direction opposite 
to the hard gluon due to the negative interference effects. In the limit when 
the hard gluon becomes softer interference is still destructive and one expects 
similar suppression. 

The lowest order result for the normalized correlation reads in the soft 
gluon limit 

C(,,^)^l + j^ ;°^^ , (102) 

26f coshr/ — cos</3 

48 



• ZEUS 95 0=60 



• ZEUS 95 G=60° 






0.5 


1 

0.5 


1 

0.5 


1 

0.5 




10<Q^<20 GeV^ 



I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ 



qSiGO GeV^ 








0.25 0.5 0.75 1 



w 



0.6 

0.4 

0.2 



-0.2 
0.6 
0.4 
0.2 


-0.2 
0.6 
0.4 
0.2 


-0.2 
0.6 
0.4 
0.2 


-0.2 



-20<Q^<100GeV^ 



10<Q^<20 GeV^ 




^'t^iiii 



■Q^lOOGeV^ 




Q^>2000 GeV^- 




0.25 0.5 0.75 1 

8 



Figure 11: ZEUS results for the scaling functions Y{e) = -/3In(f/fe)/6 (left) and Y{e) = 

2(3\nr/b (right) for various jet energies E =< Q > /2;b = 2/3-y^ln (Esin0/A). Sohd hnes 

givejjje scaling QCD predictions, namely cj(e,2) — 2, Eq. ( Pq ) and uj{e,2) — 2\/]~^l£, Eq. 

(BTj). Dotted line shows the constant as version of the asymptotic forms, Eq. (101). 



49 



1.4 

1.2 

1 

0,8 

0.6 

0.4 

0.2 



-0,2 



^8=91 GeV Dato, corrected 
Vs«91 GeV 
- JETSET MC 

Generator level 




o 



Doshed line 

Asymptotic curve 



0.4 

e = \r) 



0,6 0,8 1 

„)/fn(P0/A) 




Asymptotic Prediction 



• V9= 1 83 GeV 
dota corrected 

A Vs=183GeV 
Jetset MC 



e = lln(0/i>,2)/ln(P0/A) 



Figure 12: The same scaling functions as in FigJlll as measured by the DELPHI Collabora- 
tion for two highest jet energies. Results of the Monte Carlo simulations are also shown. 



with T] and tp being the rapidity and aziniuthal angle differences of two soft 
gluons. 

Theoretical improvements of this formula involve: i) corrections to the 
multiplicity flow, ii) corrections from the hard matrix element, and iii) resum- 
ing 0{as In ^gg) terms similar to z9i2 singularities considered in the previous 
section. Moreover, the sensitivity to the jet axis was elimina.t ed hv introduc- 
ing a concept of the energy-multiplicity-multiplicity correlations. c_L3 Additional 
weighting by the energy effectively defines the jet axis as the direction of a fast 
particle. Such a definition is simple to implement experimentally and retains 
all features of the QCD evolution. 

In Fig. |l^ theoretical calculations are compared with the MC results. In- 
deed the back-to-back configuration (if ~ 180°) is suppressed. The effect of 
non leading corrections is rather dramatic especially for the aligned gluons. 
At <y9 = TT the non leading corrections almost double the leading order result 
bringing theoretical predictions to within 10% of the MC simulation s. Mea- 



surements by OPAL and ALEPH have confirmed the above behavior .liz^LLz^ 
Experiment agrees with Monte Carlo up to 5 — 10% in the whole range of 
If which provides satisfactory support for higher order calculations and for 
the LPHD. In particular, experimental numbers for the back-to-back (in the 
azimuth) configuration read 

C(7r) = 0.787 ±0.002 ±0.004 (OPAL), 
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Figure 13: Azimuthal correlations of partons in a jet near central rapidity rj, at leading 
and next-to-leading orders for different A and eemparison with the Monte Carlo (HERWIG) 



result 
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Cin) 



0.783 ± 0.001 ± 0.016 (ALEPH), 



(103) 



where the first error is statistical and the second is systematic. This is to be 
compared with the MC result Ca/c = 0.8 and analytical improved Cnl = 0.93. 

5.4 The Critical Angle and the Local Parent Child Correspondence 

The success of perturbative QCD in describing hard processes only enhances 
our curiosity and the need to understand soft phenomena as well. It is an ex- 
perimental fact that many properties of hadronic final states are well described 
by equivalent partonic observables also at the scales close to, or even beyond, 
the standard range of applicability of perturbative QCD. This is the founda- 
tion of the Local Parton Hadron Duality hypothesis, which still inspires many 
theoreticians. Is it possible, that some properties of the soft, confining inter- 
actions are already hidden in the perturbative cascade? The only example of 
such a mechanism is provided by the pre-confinement scenario in early days of 
perturbative QCD .Ell In this section we would like to discuss another intriguing 
property of the perturbative cascade. Namely there exist a small critical angle 
which is still within the range of perturbative calculations, and below which 
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the distribut ion o f next generation of partons is equal to the distribution of 



1 90 
their parents Jifill This precise coincidence of the distributions of partons from 



different generations, which we call the Local Parent Child Correspondence, 
has some interesting consequences. 

Consider again the distribution of the relative angle 13'^' (^12, J', 6) dis- 
cussed earlier. In the double logarithmic approximation this correlation func- 
tion obey s th e integral equation which can be derived from the evolution equa- 
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tion dH) 

(104) 
The inhomogenous term dP'' is constructed from the product of 1-particle 
angular distributions and found as (i^^^(a;) = D^^\x)N{x), x = P-di2/A and 
with the multiplicity N in the cone of half-angle ^12 around the primary parton 
axis. In exponential accuracy d"' ~ N 



The solution of Eq. (104) can be written aaii-!. 



The resolvent R [^, ■;p- , ^^^^ 1 is nothing but the momentum distribution of 
the parent partons K in the jet (P, Q) with the condition, that their virtuality 
is bound from below by Kdi2 and not by the elementary cut-off Qo- This is 
equivalent to restricting their emission angle to be greater than 1912. 

It turns out that this distribution is non analytic at some small value of the 
relative ang le ??i2 = i?c- To see this consider the dependence of the integrand 



of Eq. ( 105 ) on the parent parton momentum K. The interplay between the 
decrease of the parent density R and the increase of the children densities d^^-* 
with K produces a sharp maximum which is however '(9 12 dependent. The 
important element of the running a^ analysis is that the position of this 
maximum is independent of Qq, since both distributions do not depend on 
Qq. Now, for small relative angles the above maximum shifts below the lower 
limit Qo/'di2 of the momentum integration. At that point the result changes 
non analytically (albeit very gently — analogously to the second order phase 
transitions in statistical systems). Namely, for ?9i2 bigger than -de, p^^-* is given 
( to the exponential accuracy) by the Qq independent value of the integrand 
at the maximum. For 1912 below 'dc, however, the result is given by the value 
of the integrand at the lower bound and obviously depends on Qa. Moreover, 
since d2{l) = 1 in our exponential approximation, the density of pairs below 
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"dc is given entirely by the density of parent partons at the minimal momentum 



,..(,,,). ^(_^_,:^). (106) 



The value of the critical angle derived from these considerations readsli£ 

^ _Qo/pey/^ , _go/ in(P9/A) y/^ 

for the constant and running as respectively. The transverse momentum 
associated with this angle grows weakly with the jet momentum P hence the 
perturbative description applies also in the new region Qo/ P < ??i2 < t^c- Even 
though some features of this result (e.g. the emergence of a sharp singularity) 
would be modified in more refined approximations, the existence of the two 
angular regions with different properties is genuine. In the large angle region, 
the distribution scales in e (see previous section) and is independent of Qo- 
On the other hand, below dc the distribution depends on Qq, which can be 
interpreted as the sensitivity to the soft interactions where the hadronization 
and confinement start becoming important. 

The numerical value of the critical angle is rather small at high energies, 
dc < 1° at LEP energies. It is interesting to note that the scaling in e is indeed 
best satisfied at large angles (^12 ^ dc), see Fig. n3, which is consistent with 
this analysis. 

Sensitivity to the cut-off Qq at small angles leads to an interesting predic- 
tion for the particle mass dependence of the produced pair. Assuming that the 
particle mass acts as an effective transverse momentum cutoff, we expect that 
the spectrum of pairs with large relative angles is mass independent, while for 
angles smaller than dc the density depends on the mass and is suppressed for 
higher masses. 

6 On Inclusive Properties of Heavy Quark Jets 

Until now we have discussed mainly the properties of jets in e^e~ collisions 
obtained by analyzing the full hadron data sample with contributions from 
all flavors. In the last few years intensive experimental studies have been 
performed with the heavy-quarks Q{c,b). Experiments at LEP and SLC have 
produced a wealth of new interesting results on the profiles of jets initiated 
by heavy quarks. The accuracy of measurements started to be comparable 
with that in the qq events. This is mainly related to the steady improvement 
in the heavy-quark tagging efficiencies. Further progress is expected from the 
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measurements at hadronic colliders and a future linear e+e^ collider. The 
principal physics issues of these studies are related not only to testing the 
fundamental aspects of QCD, but also to their large potential importance for 
measurements of heavy-particle properties: lifetimes, spatial oscillations of 
flavor, searching for CP-violating effects in their decays etc. 

A good understanding of 6-initiated jets is of primary interest for analysis 
of the final-state structure in ti production processes. A detailed knowledge of 
the 6-jet profile is also essential for the Higgs search strategy. 

The physics of heavy-quarks has been traditionally considered as one of 
the best testing grounds for QCD. Note that since the mass of 6-quark is 
much larger than the QCD scale A, the 6-quark fragmentation proves to be 
an especially useful tool for studying the perturbative scenario. In this sec- 
tion we will restrict ourselves to the discussion of some specific features of 
heavy-quark-initiated events which are related to our previous considerations 
of multiparticle production and LPHD. 

6.1 QCD Bremsstrahlung Accompanying Heavy Quark Production 

When examining the perturbatively based structure of the heavy-quark jets, 
one has to understand first how the development of the parton casca des gen- 
erated by Q depends on the quark mass M. It was Hpmrmt^tratprlP'-H^M'l^^^l^-^o 
that the difference in many properties of hadronic jets initiated by heavy quarks 
(excluding the products of the weak decay of the heavy quark iteself), from 
that of light quarks, originates from the restriction of the phase space avail- 
able to gluon radiation associated with the kinematic effects of the heavy-quark 
mass. 

As well known, the radiation pattern of the primary soft gluon with energy 
u) from a massive relativistic quark with energy Eq 3> M and small emission 
angle 9 ^ 1 is given by 

, as „ Q^dQ^ du; 

where 

M 
90 - — . (109) 



From Eq. (108) one concludes that the large double-logarithmic contribution 



comes only from the region of relatively large radiation angles 

9 > 9o, (110) 
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where emission becomes insensitive to M and appears to be identical to that 
for the hght g-jet 

as ^ d&^ duj 

For the region 



d(TQ^Q+g = -^Cf -^^ -7. (Ill) 



e < Go (112) 

the angular integration is no longer logarithmic. 



As it follows from (|108|) the forward gluon radiation is suppressed .1 ^1 '^ 
This phenomenon is characteristic for bremsstrahlung off a massive particle. It 
reflects the conservation of the projection of the total angular momentum on 
the massive fermion momentum: the soft radiation does not change the quark 
helicity and the forward emission is forbidden (in analogy to the well-known 
"0 — transition" phenomenon). As it is discussed below, after taking into 
account parton cascades generated by a primary gluon, the relative yield of final 
particles due to this part of phase space is 0{^/as). Since the main logarithmic 
contribution is absent and the differential cross section ( |108| ) vanishes in the 
forward direction, it is naturalMrl to call this region the "dead cone" — a 
relatively depopulated cone around the quark direction with an opening angle 
6^60. 

Notice that the standard expression for the soft bremsstrahlung spectrum 

2CFas v^ sin^ 6 dtj 
d<7Q^Q+g = 7; 5 5777 — dcosB (113) 

■K [1 ~ v' COS^ B)^ UJ 



with 




(114) 



describes adequately the non-logarithmic region of finite emission angles B ~ 1. 
This part of the radiation phase space contributes as ©(y^os) too. Correspond- 
ing corrections will be taken into full account later when we derive formulae 
for particle multiplicities outside the scope of the double logarithmic approxi- 
mation. 



Recall that the angular pattern in (113) is exactly the same as in t he Q ED 



1 9f 
case of classical photon bremsstrahlung off the massive electric charge .Iir2; 

As we mentioned above, the structure of the primary gluon radiation at 



large angles (110) appears to be identical to that for the light q-jet. The same 
holds true for the internal structure of secondary gluon subjects (angular- 
ordered cascades). 
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6.2 Specific Features of Events Containing Heavy Quarks 

As was discussed in detail in the previous sections, the dominant role of the 
perturbative dynamics has been very successfully tested in studies of light- 
hadron distributions in QCD jets. Within the LPHD approach the quark 
mass-induced restriction on the gluon radiation leads to a number of important 
consequences for the Q-quark jets. 

(i) Explicit visualization of the dead cone 

From the measurements of the "companion" particle distributions inside fixed 
cones around the Q direction one can learn about the expected depletion of 
the forward-particle production for angles 9 < J^ . Note that a study of the 
forward hadroproduction in Q-jets which is free from "trivial" perturbative 
bremsstrahlung effects, may help in the investigation of the subtle features of 
some non-trivial confinement scenarios. 

The first comparison of the primary particle angular distribution in t he b 



and light-quark jets in e+e" annihilation has been reported by DELPHI.LL£l 
In the data analysis special care has been taken in order to exclude the decay 
products of the 6-hadrons from consideration. Preliminary results look quite 
encouraging. However further detailed studies are needed for establishing the 
perturbative-QCD origin of the observed depopulation of the forward-particle 
production in the 6-quark initiated jets. 

(ii) Leading particle effect in heavy-quark fragmentation 

Suppression of small angles restricts emission of hard gluons with relatively 
small transverse momenta k±, where radiation is normally most intensive as 
being proportional to as{k'j^). This results in a decrease in the energy frac- 
tion tha tj th(; hnavy -qiifi,rk Q is sharing with the secondary bremsstrahlung 



particles li^y.^1^ Thus, the inclusive distribution of Q is expected to peak 
near xq ~ 1, the so-called leading partic le ef fect, which first was predicted 
within the framework n f thp parton modelliHrl and quantified as the so-called 



1 "^n 

Peterson fragmentation.li21J Recall that at the hadronization stage the heavy- 



quark loses only a momentum fraction of order A/M when forming a heavy- 



light hadronJizj The leading heavy particle phenomenon has been cle arlv ob 



L31132 



served experimentally in the bb and cc events in e'^e" annihilation 

As expected, at the Z° the 6-quark inclusive spectra peak very close to one, 

around xq ~ 0.8 — 0.9. 

The record accuracy in determinatio n of the &-quark spectrum has been 



achieved in the recent SLD measurement .lL2iJ As a result, this allows to dis 



criminate between current parton-shower plus hadronization models. 
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As discussed in detail,Ll-£j within the LPHD scenario the effects of pertur- 
bative gluon radiation allow to reproduce the shape of the heavy-quark inclu- 
sive spectrum provided the perturbative prediction is extrapolated smoothly 
down to the region of small gluon transverse momenta within the universal low- 
scale as approach. As well known such approach proves to be quite successful 
in the description of the jet event shape, for a discussion, see for example, 
Ref. 13 and the article by V. Braun and M. Beneke in this book. 

The Sudakov form- factor suppression of the quasi-elastic region at xq — > 1 
results in the distribution which is qualitat ively similar to the parton-model- 

h on] 

motivated Peterson fragmentation function.t2!:ll Note that a,1l- orde r resumma- 
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tion of the coUinear and Sudakov logs have been considered alsoJiji3 It is worth- 
while to mention, that as far as the heavy-quark is concerned, based on LPHD 
one may expect that a purely perturbative treatment is dual to the sum over all 
possible hadronic excitations. Therefore, without involving any phenomeno- 
logical fragmentation function at the hadronization stage, one could attempt to 
describe the energy fraction distribution averaged over heavy-flavored hadron 
states, the mixture that often appears experimentally. Such purely perturba- 
tive approach is at least free from the problem of "double counting" which one 
may face when trying to combine the effects of perturbative and hadronization 
stages. 

(iii) Average multiplicity of events containing heavy quarks 

As we discussed above, a forward suppression of soft-gluon radiation off an 
energetic massive quark Q induces essential differences in the structure of the 
accompanying radiation in light- and heavy-quark-initiated jets. According to 

the LPHD concept, this, in particular, should lead tp-LDnjcapanding difference 

1 9^ 1 '-{4 ^Kl 
in companion multiplicity of radiated light hadrons J££2^K '^'W ^ 

It is a direct consequence of the perturbative approach that the difference 
of companion average multiplicities of hadrons, AA^g^ from equal energy (hard- 
ness) heavy- and light-quark jets should be energy-independent (up to power 
correction terms 0{A'P /W'^). This constant is different for c and b quarks and 
depends on the type of light hadron under study (e.g. all charged, tt*^, etc.). 
This is in marked contrast with t he p rediction of the models based on the idea 



of reduction of the energy scale,L22l Nqq{W) ~ Nqq{{l — {xq))W), so that 
the difference of q- and Q-induced multiplicities grows with W proportional to 
N{W). We shall elucidate this bright prediction of the perturbative scenario^ 

'The discussion of companion multiplicity presented below was borrowed from the pa- 
per by Yu.L. Dokshitzer, S.I. Troyan and one of the authors (VAK) which has never been 
published. 
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Recall, that companion multiplicity is an infrared-sensitive quantity dom- 
inated by the emission of relatively soft gluons. In this domain the only dif- 
ference between heavy- and massless-quark cases comes from the suppression 
of soft bremsstrahlung in the forward direction 81+ < 80 = M/Eq, where 
81+ is the angle between a primary gluon [gi) and the quark (+). In the 
region of large gluon radiation angles 81+ 3> 80, the finite-mass effects are 
power-suppressed and do not affect the picture of strictly angular-ordered evo- 
lution of gi as a secondary jet. It is straightforward to verify that within the 
MLLA accuracy the finite mass induces only a small integral correction to the 
angular-ordered prescription for the next-generation gluon 52. As a result, in 
the large-angle region the internal structure of secondary-gluon jets is identi- 
cal to that for the light-g case, and the emission angle 81+ should be taken 
as an evolution parameter to restrict the subsequent cascading. As usual, the 
smaller the radiation angle 81+, the less populated with offspring partons the 
gluon subjet gi is. 

The situation changes, however, when 81+ becomes smaller than 80. The 
dead cone region gives a sizeable (though nonlcading) contribution, of order 
^/os, and should be taken into account within the MLLA. Here the opening 
angle of the jet 171 "freezes" at the value 80 and no longer decreases with 
81+ -^ 0. The reason for this is rather simple. Normally, in the "disordered" 
angular kinematics 821 > 81+ the destructive interference between emission 
of a soft gluon 52 by quark and gi cancels the independent radiation gi -^ g2- 

Meantime, in the massive-quark case the interference contribution enters 
the game only when 82+ > 80, so that the cancellation of the independent 
1^2 term inside the 80-cone does not occur. In physical terms what happens 
is the loss of coherence between -|- and 1 as emitters of the soft gluon 2 due 

(2) —1 

to accumulated longitudinal separation Azi_|_ > Am « 0^2 between massless 

and massive charges (wi = l,f+ ~ 1 — 89/2 < 1). Indeed, during the for- 

(2) 
mat ion time of g2itf ^ 

longitudinal direction by 



niation time of g2,tf ^ 1/0^2821 the quark and the gluon 1 separate in the 



Azi+^if l^+-cos8i+| ^ X\f ■{Ql+ + el)/elv (115) 

It is the last factor which determines whether interference is essential or not. 
When this ratio is larger than unity, Q and gi are separated enough for 52 to be 
able to resolve them as two individual classical charges. In these circumstances 
gluon gi acts as an independent source of the next-generation bremsstrahlung 
quanta. Otherwise, no additional particles triggered by gi emerge on top of 
the yield determined by the quark charge (which equals total color charge of 
the Qgi system). In the massless limit (80 — 0) this reproduces the standard 
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angular-ordered picture, O21 < Qi+. In the massive case one concludes from 



(115) that the upper limit of the relative gluon angle remains finite when 01+ 
falls inside the dead cone: 621 < 6o for arbitrarily small 9i+ ^ 6o- 

Thus, the proper evolution parameter for subsequent parton cascading of 
the primary gluon gi (generalization of an "opening angle" of jet 1) may be 
chosen as 

el+ ^ el+ + el (iie) 

Another comment is in order concerning generalization of the argument of the 
run ning coupling to the massive-quark case. Here again the substitution similar 
to ( |116| ) is applicable. Namely, for the effective coupling that determines the 
probability of the emission Q —^ Q + gi one has to use as an argumentlHil 



kj_ = uj 



2 sin 



61- 






Lole 



!+■ 



(117) 



Within MLLA the expression for the multiplicity of light particles ac- 
companying the production of a heavy-quark pair, Nqq, can be obtained by 
convoluting the probability dw of a single gluon bremsstrahlung off a heavy- 
quark Q with the parton multiplicity initiated by the gluon subjet with the 
hardness parameter k±. 

The differential first-order cross section for e+e^ -^ QQ a integrated over 



the angles of the QQg system was f irst presented by Ioffdi22l (QED first-order 



results can be found elsewherdi^il) . For the case of vector QQ-production 
current the differential over the gluon energy fraction z 



E 



= -i, 2E = W = 



gluon emission probability can be written aslizSlir^ 



dw\ 



Cpcts dz drj 



VT^ 



2(1 



Here 



1 > ry 



1 



[i^cos^Qc > 






Am^ 



+ z' 



(118) 



(119) 



% 



1 



m = MjW < 1, 



C^y = (3-w2)/2 = 1 -t- 2m2 (120) 

with /? the quark velocity and Oc the polar gluon angle in the QQ rest frame, 

Arr? 



a' = 0\z) 



1 



l-z 
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< v' 



1 



4m > z. 



(121) 



The first term in curly brackets contains the main (double logarithmic) con- 
tribution and corresponds to the universal soft gluon bremssirahlung. As it 
follows from the celebrated Low soft bremsstrahlung theorem,li^ extended to 



the fermion case by Burnett and KroUJiiiil both z ^ and z° parts of the radi 



ation densit y hav e classical origin and are therefore process (and quark spin) 
independent .Ii22l This first term explicitly exhibits the dead cone behaviour: 
forward "soft" radiation vanishes, r\ ^^ r]Q when sinSc — » 0. The second term 
which is proportional to z^ (hard gluons) depends, in principle, on the QQ 
production mechanism. Namely, both — ^ subtraction term a nd the C- factor 

may be different for the different production currents. lijiJiizj This indicates 
that already at the level of as corrections the average companion multiplicity 
acquires a process-dependent contribution from the three-jet ensembles and, 
thus, cannot be treated any more as an intrinsic characteristic of the QQ sys- 
tem. 

Neglecting now the corrections of relative order 0.3 and rri} the average 
companion multiplicity 7Vnn(M^) can be given by 



Nqq{'^) 



■K 



dz 

z 



Vn 



dr] 1 + (1 - zf 



1^2(1 



z)\ (asik^) Ncik^)) 



(122) 



with fc^ = {zW/2)'^ri. Note that the second term in the curly brackets is non- 
logarithmic in 7] and generates next-to-leading correction SiN ^ ^J as{M)N {M) , 
see below. 

The main contribution to Nqq comes from the DL phase-space region. 
Meanwhile, the expression (122) keeps track of significant SL effects as well, 
provided that the multiplicity factor Nq is calculated with the MLLA accuracy. 

Introducing a convenient variable 






one can rewrite Eq. ( |122| ) as 

Nqq{W) 



AP 



df r 
.2 



K^ 



1 



E' 






2 sin 



Oh 






(123) 



dk I 



Qo 



P^iz)^^NGik^), (124) 
Zn 



with k± = ZK. Here P^ stands for the standard DGLAP kernel and Qo de- 
notes a transverse momentum cut-off as used in the previous sections. In the 
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massless limit M < Qq the contribution of the A'P / n^ term vanishes and (124) 
reproduces the known result for the q-']ct multiplicity. The integration can be 
performed with the use of the relation corresponding to the evolution equation 
for light-quark jet multiplicity: 



with 



N,,{Q) = 27Vq ( I ) , (126) 



see Section 3. Then Eq. (124) takes the form 



Nqq{W) = N,,{W) - 2N, (^^^ - 2N'^ (^^) + OiasN,;QlN'^). (127) 

Notice that the factor 2 in the argument of A^g generates ^/asNq{M) correction 
and is under control in the present analysis whereas it could be omitted in the 
N' t erms as producing asNq{M) terms which we neglected systematically in 
(124). Within this accuracy the term N' ~ ^JasNq{M) can be embodied into 



the multiplicity factor by shifting its argument, namely 

„, (f ) . K (f ) ^ ». (m^) . (-. 

Finally we arrive at the formula expressing the companion multiplicity in 
e+e^ -^ QQ in terms of that of e+e^ -^ qq (assuming M 3> A) 

Nqq{W) = Nqq{W) ~ Nqqi^M) [1 + 0{as (M))]. (129) 

The total particle multiplicity in QQ events then reads as 



»+ 



N^ - -W(^) ^ Nqq{W) + 7l^^ (130) 

where n'^ stands for the constant decay multiplicity of th e heavy-quarks in'^ — 
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11.0 ± 0.2 for 6-quarks, n^'^ = 5.2 ± 0.3 for c-quarkti!j2a) . 

Some comments are in order here. The main consequence of ( |l30| ) is that 
the difference between particle yields from q- and Q-jets at fixed energy W 



depends on the heavy-quark mass and remains VF-independent.B_jilJL_2zl The 



question arises, to what extent this difference can be quantitatively predicted? 
An important issue concerns non-leading corrections '^ OioLgiy^y) to the main 
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term. From the first sight one may expect that instead of (129) the correct 
formula would look like 

Nqq{W) ^ {N,,{W)}^'^'^^[l + 0{as{W))] 

- {7V,,(V^Af )}^i^LA [^ ^ 0{as{M))], (131) 

where {jV}'^^'^^ denotes the MLLA multiplicities accounting for the ^/as + as 
effects in the anomalous dimension and 1 + y/as terms in the normalization 
(coefficient function). Meantime, because of the steep growth of the multi- 
plicity (faster than any power of InW) the neglected, order ^/as(W)N{W) 



terms would dominate over the finite subtraction term in (131) with increasing 
energy. Thus, the very possibility to discriminate between the Q- and q-quark 
jets within the present theoretical accuracy may be endangered. However, a 



close inspection of the subleading corrections to Eq. (129) (proportional to 
Nqq{W)) shows that all of them are independent of the heavy-quark mass M. 
Thus, the first corrections of order as{W)Ngg{W) arise either from further 
improvement of the description of anomalous dimension Aj^as) ~ Og deter- 
mining intra-jet cascades, or from 0{as{W)) terms in the coefficient function 
due to 

• 3-jet configuration quark + antiquark + hard gluon at large angles, 

• the so-called "dipole correction" to the angular-orderin g sc heme quark + 
antiquark -f two soft gluons with large emission anglesf^ 

which are insensitive to the Oq value with the power accuracy ~ 0q ^ 1. 



Therefore, replacing the approximate MLLA multiplicity factors in (131) by 



the real observable multiplicities, one arrives at (129) which makes it possible 
to establish a phenomenological relation between light and heavy-quark jets 
with relative accuracy ^Jas{M) M'^/W'^. 

Thus, the difference between particle yields from o- and Q-jets at fixed 



annihilation energy W proves to be W^-independent P^lr^^l "^ 

(5q, = N'^^^'-'^'^iW) - N^^^'^o^W) ^ const (W), 



^bc 



f>'T' P ^Ifl) /-r-r-r\ ■»-r/:.+ 



N" " ^""{W) - N" " ^""{W) = const (M/). (132) 



Let us emphasize that it is the QCD coherence which plays a fundamental 
role in the derivation of this result. Due to this the gluon bremsstrahlung 
off massive and massless quarks differ only at parametrically small angles 
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Figure 14: The difference &9 between the average charged multipUcities of b- and hght-quark 



events 
of (SP^ 
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bq 



The perturbative expectation is chosen as the average of the experimental values 
jj up to the Z^ . Also shown is the prediction of the model based on the reduction of 
the energy scale (independently of the quark flavor). 



< 8o = M/E where, due to the angular-ordering, cascading effects are 
majorated by the N'{M) factor. 

The results of experiments on multiplicities in bb and cc events of e^e~ 



annihilation 



141 



u| 



show that to the available accuracy the differences dhq and 
5cq are fairly independent of W . This is in marked contrast to the steeply 
rising total multiplicity and confirms the MLLA expectations. As i llustrated 
in Fig. n4 the data on 5bq are cl early inconsistent with a mode] 
the reduction of the energy scale. 



135 



based on 
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JTc/l 



We turn now to the absolute values of the charged multiplicity difference 



Neglecting in Eq. (12£) t he s ubleading 0{as) corrections, the following 



1 ^4 

MLLA expectation was foundli2zl by inserting the appropriate experimental 



numbers into the r.h.s. of Eqs. (129)-(130) 



51^ = 5.5 ± 0.8 



(133) 



This exceeds the experimental values in Fig. O, showing the essential role of 



the next-to- MLLA (order as{M) ■ N{M)) terms. An attempt was madellirl to 



improve Eq. (|129|). However, we have to mention here that the very pict ure of 



1 4*^ 
accompanying multiplicity being induced by a single cascading gluorii-i2 



is not 
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applicable at the level of subleading effects. Therefore, a self-consistent reliable 
theoretical improvement of the MLLA predictions for the absolute values of 
Nqq remains to be achieved. 

Further detailed experimental results could provide stringent tests of the 
perturbative predictions. It will be very interesting to check whether the differ- 
ences 5^ of mean multiplicities for the identified particles remain energy inde- 
pendent. Results from the Z^ are available as reference values. In particular, 
t he D ELPHI dataliz3 suggest S^ ~ and therefore, using again Eqs. (129)- 
(|130|), one may expect the relation 



N^o^ ^""/Ny -*«« ~ 1. (134) 

The energy behaviour of 5q should be watched closely at a future linear 
collider. 

(iv) Spectra of light particles from the heavy- quark jets 

The bremsstrahlung-suppression effect on the energy spectra of light particles 
accompanying QQ production can be examined in a very similar way to the 
average multiplicity considered above. The resulting companion particle distri- 
bution DQ{x,ln{E/A)) is predicted to be depopulated in the hard momentum 



region compared with the case of the light-quark I?^(a;,ln_B/A).lir£lcy It is 



instructive to first gain insight by considering the double- logarithmic result. 
Due to the dead-cone phenomenon the difference between Dq and D^ is con- 
nected with the gluon radiation at angles Q < Qo- Radiation in the restricted 
angular cone is similar to the 9-jet production process with the characteristic 
hardness M ^ A. This process induces a gluon radiation pattern with the 
corresponding cascading which leads to a certain final system of light hadrons. 
Now, boosting this jet with the Lorentz-factor 7 — E/M along the direction 
of its momentum one arrives at an exact image of the dead cone: an ensemble 
of energetic particles (with energies up to E) concentrated inside the cone 80 . 
Thus, with the double logarithmiq_acciixacy a simple formula for the Q-jet 
particle distribution can be written 



iwm 



^h I ^ ,„ E\ _ . ( E\ ,^ f , M 



D'q [^M jj = D- i^x,\n -j D- (^x,ln -j . (135) 

Here the structure of the heavy-quark jet is expressed in terms of the particle 
distribution generated by the light-quark at different energy scales: Q" and 



AP. In Eq. (135) x is the light-hadron energy fraction, x = Eh/E > A/M. 
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Therefore the depopulation occurs for quite energetic particles only, whereas 
in the soft domain 

K<E,,<E^^ (136) 

the light-hadron spectra in the Q- and light-quark jets should be identical. 

In order to obtain the MLLA result one has to take into account, first of 
all, the energy loss by heavy-quark at the first steps of the evolution. Using the 
evolution equations the expression for the Q-jet companion particle spectrum 
can be approximately given by 

where D'^{x,\wW/A) is the standard spectrum of particles in a g-jet and {y) 
being the averaged scaled energy fraction of the Q. This problem needs further 
detailed studies. 

Measurements of t he in clusive charged particle distributions in b- and uds- 
quark events at the Z'lizB clearly show that for 6-quarks the momentum spec- 
trum is significantly softer than for uds quarks. Unfortunately so far no com- 
parison with the analytical QCD expectations has been performed. 

7 Color-Related Phenomena in Multi-Jet Events 

7. 1 Radiophysics of Hadronic Flows 

Coherence phenomenon is an intrinsic property of QCD (and in fact of any 
gauge theory). Its observation is important in the study of strong interactions 
and in the search for deviations from the Standard Model. As we have already 
discussed, within the LPHD scenario, an intimate connection is expected be- 
tween the multi-jet event structure and the underlying color dynamics at small 
distances. The detailed features of the parton-shower system, such as the flow 
of color, governs the distribution of color-singlet hadrons in the final state (see 
e.g. Refs. 15,146). 

Color coherence effects in hadron multiplicity flows in the inter-jet regions 
have been very well established from the early 1980's in e+e^ annih ilati on. 



see Refs. 147,148,106,149,150, in what has been termed the "string"E^ 

E4I 
"drag"0 effect. As discussed below the existing data are clearly in favour of 

QCD coherence, in particular, particle production in the region between the 

quark and antiquark jets in e^e~ — > qqg events is suppressed. In pQCD such 

effect s ar ise from interference between the soft gluons radiated from the q, q 



and (7JHSIL2I The PETRA/PEP data first convincingly demonstrated that the 
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wide-angle particles do not belong to any particular jet but have emission prop- 
erties dependent on the overall jet ensemble. The inter-jet-coherence phenom- 
ena were then successfully studied at LEP, TRISTAN and TEVATRON (also 
discussed in recent reviews! '11^^ rr^^). The experiments have nicely demon- 



strated the connection between color and hadronic flows. 

Surely, it is entirely unremarkable that the quantum mechanical interfer- 
ence effects are observed in QCD. Of real importance is that the experiment 
proves that these effects survive the hadronization phase. 

The inter-jet coherence deals with the angular structure of particle flow 
when three or more hard partons are involved. The hadron distribution proves 
to depend upon the geometry and color topology of the hard-parton skeleton. 
The clear observation of inter-jet-interference effects gives another strong ev- 
idence in favor of the LPHD concept. The collective nature of multiparticle 
production reveals itself here via the QCD wave properties of the particle flows. 

The detailed experimental studies of the color-related effects are of par- 
ticular interest for better understanding of the dynamics of hadroproduction 
in the multi-jet events. For instance, under special conditions some sub- 



tle interfemn 
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c.p e ffects breaking the probabilistic picture, may even become 
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dominant .ti^li23 We recall that QCD radiophysics predicts both attractive and 
repulsive forces between the active partons in the event. Normally the repul- 
sion effects are small, but in the case of color-suppressed 0{1/Nl^) phenomena 
they may play a leading role. It should be noted that usually the inter-jet 
drag effects are viewed only on a completely inclusive basis, when all the con 



stituents of the multi-element color antenna are simultaneously active E^'IL^Z 
A challenging possibility to operate within the perturbative scenario with 
the complete collective picture of an individual event (at least at very high 
energies) was first discussed in Ref. 155. The topologometry on the event- 
by-event basis could turn out to be more inform ative than the results of mea- 
surements averaged over the events.li23lL±lllL22l Note, that there is an essential 
difference between the perturbative radiophysics and the parton-shower Monte 
Carlo models. The latter not only allow but even require a completely exclusive 
probabilistic description. Normally (such as in the case of e^e~ — » qqg) the 
two pictures work in a quite peaceful coexistence; the difference only becomes 
drastic when one deals with the small color-suppressed effects. 

Note that in spite of theoretical uncertainties, the phenomenological suc- 
cess of the Monte Carlo models which include interfering gluons indicates that 
the color-suppressed interference terms do not induce a very large effect. 

Let us emphasize that the relative smallness of the non-classical effects by 
no means diminishes their importance. This consequence of QCD radiophysics 
is a serious warning against the traditional ideas of independently evolving par- 
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tonic subsystems. So far (despite the persistent pressure from the theorists) 
no clear evidence has been found experimentaUy in favour of the non-classical 
color-suppressed effects in jets, and the peaceful coexistence between the per- 
turbative inter-jet coherence and color-topology-based fragmentation models 
remains unbroken. However, these days the color-suppressed interference ef- 
fects attract increased attention. This is partly boosted by the findings that the 
QCD interference (interconnection) between the W~^ and W~ hadroni c decays 
could affect the precise W mass reconstruction at LEP-2, for example. 
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QCD intercnnnectio ri ma y affect also the top quark studies, in particular, its 



mass reconstruction. 
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It is worthwhile to mention that the relative smallness of the color-suppressed 
interference effects is supported by recent s earc hes by OPAL of the so-called 
reconnection effects in hadronic Z^ decays.Lt£i] Finally, we recall that color- 
related collective effects could become a phenomenon of large potential value 



as a new tool helping to i ^istjpp;ii^|sh 
tional QCD backgrounds 



128 146 15 



ho now phy sics signals from the conven- 
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In the next subsections we briefly survey the basic ideas related to QCD 
collective effects in multi-jet events and then we describe (wherever possible) 
the latest data on QCD radiophysics. 



7.2 String/Drag Effect in qqg Events 

The first (and still best) example of the inter-jet color-related phenomena is the 
string/drag effect in e"'"e^ -^ qqg. Since the new results based on the refined 
analysis of qqg events continue to pour out from the LEP groups it seems useful 
to recall the main ideasO of the pQCD explanation of this bright coherence 
phenomenon. We consider first the angular distribution of particle flows at 
large angles to the jets in e+e^ -^ qqg. The more general case including the 
inside jet particle flow will be discussed below. 

Let all the angles 8^^ between jets and the jet energies Ei be large (i = 
{+ - 1} = {qqg}) : 9+- - G+i - e_i - 1, £;i - £;+ - £;_ - £; - W/i. As 
was discussed above, within the perturbative picture the angular distribution 
of soft inter-jet hadrons carries information about the coherent giuon radiation 
off the color antenna formed by three emitters {q,q and g). The wide-angle 
distribution of a secondary soft gluon 172 displayed in Fig. n_q can be written as 



SndNgg 

d^n2 



— WMn2)N'^iYrr.) = 



(1^) + (i^) - ^ (+^) 



(138) 
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Figure 15: Kinematics of inter-jet radiation in three-jet events. 



Here the "antenna" {ij) is represented as 



m 



aitti 



Oii = (1 - nifij), 



(1 - n2ni) 



(139) 



and N' (Ym) is the so-called cascading factor taking into account that a final 
soft particle is a part of cascade (see Refs. 64,15 and Section 3), Ng{Y) = 
dNg/dY. Furthermore, Ym = lnii^8m/A, where one defines the angle Qm = 
min{0+, 8_, &i} with cosOi = n2ni for i = {-I-, — , 1}. 

The radiation pattern corresponding to the case when a photon 7 is emitted 
instead of a gluon reads 



SirdNq, 
dQft2 



— W+^{n2) N'g{Y^) = -^ (+^) 7V;(y„). (140) 



The dashed line in Fig. |l^ displays the corresponding "directivity diagram" , 
which represents the particle density (140) projected onto the 557 plane: 



V{a,P) 



2Cf 



d cos 02 



cos a — cos P \ sin a 



(+-) = 2CFa+^V{a,P) 
- ~ a TT — /3 



sin/3 



1^92, /3 = 9 + - - ¥52- 

(141) 



The expression VF-^ (n^) is simply related to the particle distribution in two- 
jet events e^ e~ -^ <l{p+) + <i{P-)j Lorentz boosted from the quark cms to the 
lab system. 

Replacing 7 by gi changes the directivity diagram essentially because the 
antenna element gi now participates in the emission as well. However, this 
change leads not only to an appearance of an additional particle flow in the gi 
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Figure 16: Directivity diagram of soft gluon radiation, proj ected o nto the qq^y (dashed line) 
and qqg (solid) event planes. Particle flows of Eqs. (1411) and ( 145 ) are drawn in polar 
coordinates: = ip2,r = \n2W{ip2). Dotted circles show the constant levels of particle 

flow: W{ip2) = 1,2,4. 



direction. Integrating (|138|) over 02 one obtains (7 = O^ 



^2) 



W±i(ip2 



Nc 



i+i 



l/(a,7) + a-iT/(/3,7) 



Nh 



a+-Via,P) 



(142) 



Fig. nfl illustrates that the particle flow in the direction opposite to Hi ap- 
pears to be considerably lower than in the photon case. So, the destructive 
interference diminishes radiation in the region between the quark jets giving 
a surplus of radiation in the q — g and q — g valleys. One easily sees that the 
color-coherence phenomena strongly affect the total three-dimensional shape 
of particle flows in three-jet events, practically excluding the very possibility 
of representing it as a sum of three parton contributions. 

Let us recall here that, owing to coherence, the radiation of a soft gluon 32 
at angles larger than the characteristic angular size of each parton jet proves 
to be insensitive to the jet internal structure: 52 is emitted by a color current 
which is conserved when the jet splits. This is the reason why one may replace 
each jet by its parent parton with pf w 0. 

For illustration of the particle drag phenomenon it is instructivel2j to exam- 
ine the fully symmetric qqg events (Mercedes-type topology) where n^n^ = 
n+fii = n-fii. Let us take n2 pointing in the direction opposite to Hi, i.e. 
midway between quarks. Thus, neglecting the weak dependence N' on 0, one 
arrives at the ratio (numbers for Nc = 3) 



dNqqg/dn2 
dNqq^/dn2 



Nl-2 



2{Nl 



0.44. 



(143) 
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The corresponding ratio for the projected particle flows obtained from (141) 
and (|l4|) reads 

dN,,^/dip2 14(7V2-1) "•°"- ^'^^^ 

We emphasize that Eq. ( |13^ ) provides not only the planar picture, but the 
global three-dimensional wide-angle pattern of particle flows. 

For hadron states associated with the qqg and qqj events the above ratios 
of particle flows should remain asymptotically correct, since non-perturbative 
hadronization effects are expected to cancel at high energies. It is an interesting 
question, to what extent the non-perturbative effects are important at present 
energies for a quantitative description of the data. 

The analysis of the bremsstrahlung pattern clearly demonstrates particle 
"drag" by the gluon jet gi. If one drops the color suppressed contribution, 



the two remaining terms in Eq. (138) may be interpreted as the sum of two 
independent (1-I-) and (1—) antenna patterns, boosted from their respective 
rest frames into the overall qqg cms. The point is, that by neglecting the 
1/N^ terms, the hard gluon can be treated as a quark-antiquark pair. In this 
approximation each external quark line is uniquely connected to an external 
antiquark line of the same color, forming colorless qq antennae. In the gen- 
eral case, when calculating the resultant soft-radiation pattern, one can only 
deal with a set of such color-connected qq pairs because the interference be- 
tween gluons emitt ed frorri nnn-cn lnr- connected lines proves to be suppressed 
by powers of i/iV^ .lllislMllM _ 



Note that in the first order in as the qq-antenna pattern (140) can be 
presented in the form^ 

, ^ as{k±) dk\ 
dwqq^ = Cf TJT "J^ll (145) 

TT rZ I 

where 

y„ - \ ln^-|4, fci - '^--'^^^P-'^^ (146) 

2 (p--fc2) {P+-P-) 

are the Lorentz-invariant generalizations of the rapidity and the transverse 



momentum of (72 in the qq cms correspondingly. Eq. (145) describes the struc- 
ture of the gluon radiation in e+e^ — > qq events. Due to color coherence the 
particle multiplicity in a secondary gluon jet (72 dep ends on the transverse mo- 



mentum k I of t he gluon and not on its energyJ^'E2l The radiation a nten na 



formalismoli^ is closely related to the Lund color dipole approach.li^ In 
both schemes the natural choice for the cut-off parameter in the cascades is the 
gluon transverse momentum in the cms of the emitting antenna-dipole. This 
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implies that the soft gluons connect the active hard partonsjn, exactly the same 



way as the string in the Lund string fragmentation model |^ I which illustrates 
the connection between pQCD and the string modelEj Notice, however, that 
within the latter model there is no string piece spanned directly between the 
quark and antiquark, no particles are produced in between them, except by 
some minimum "leakage" from the other two regions. 

When a gluon is emitted in, let us say, an e'^e^ annihilation event, then 
the radiation of a second softer gluon is described by two antennas-dipoles, one 
spanned between the q and g and another between the g and q. The emission 
of a third, still softer, gluon is described in terms of three antennas-dipoles etc. 
The gluons appear to be color-ordered in such a way that the dipole is stretched 
between the color charge of one gluon and the corresponding anticolor charge 
of the next. Thus, neglecting the l/N^ terms, the QCD cascade is formulated 
as a branching process in which antennas-dipoles are successfully split into the 
smaller and smaller ones. 



7.3 Multiplicity Flows in Three- Jet Events 

We begin here by comparing particle flows in qqg and qqj events and in the 
next subsection we consider the respective total event multiplicities. 

In what follows we first do not refer to a particular three-jet-event selection 
method but focus nu th e soft radiation by the qqg system for three given 



angular directions 
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In the leading order in as the massless parton kinematics is unambiguously 
fixed as follows: 



sinGi- „ sin 61+ sin 6+ 

x+ — 2 — — ; — — — , X- — 2 — — ; — — — , Xl = 2 



E sin Qij E sin Qij E sin Oy 



X, =2, 



(147) 



with Xi = 2Ei/W being the normalized parton energies and Qij the angles 
between partons i and j{+, — = q, g; 1 = gi). We emphasize here, that, owing 
to color-coherence, the radiation of a secondary gluon 52(^2 <C Ei) at angles 
larger than the aperture of each parton jet is insensitive to the jet internal 
structure. 

Let us turn to the radiation pattern for qqj events. In the previous sub- 
section we considered the radiation at large angles to the jet directions, here 
we include the small-angle radiation in MLLA. The angular distribution of 
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particle flow can be written as 



B7rdiV„ 



d^ri2 



— N'g{Y,+ ,Y,) + — N'^{Y,^,Y,) + 2/+_ N'^{Y), (148) 



where 



Yq(q) = bl- 



and 



^q{q) 



A 



Y„+ = In 



£;,ya772 



A 



Ys_ = In 



E^^^r/2 



A 



r = ln 



E 



I+- = (+-) - 



A 
(149) 

(150) 



The factor N'^{Yi, Y) = {d/dYi)NA{Y^, Y) takes into account that the final 
registered hadron is a part of cascade. NaIYi^Y) stands for the multiplicity 
in a jet A{A = q, g) with the hardness scale Y of particles concentrated in the 
cone with an angular aperture Qi around the jet direction n^. To understand 
the meaning of the quantity NA{Yi, Y) it is helpful to represent it as 



NAiY,,Y) = 
AC = 



E 

B=q,g 



dzzDA{z,AONB{Y,), 



in(r/r. 



Yi = Fi + In z = In 



zE 

X 



.(151) 



Here NsiYi) is the multiplicity in a jet with the hardness scale Yi, initiated 
by a parton B within the cone 8i, and Da denotes the structure function for 
parton fr agm entation A ^ B. 

Eq. (151) accounts for the fact that due to the inside-jet coherence the 
radiation at small angles 8i ^ 1 is governed not by the overall color of a jet 
A, but by that of a sub-jet i?, developing inside a much narrower cone Qi. 

For the emission at large angles (a+ ^ a_ ^ 1) when all the factors N' 
are approximately the same, Eq. (|148|) reads, c.f. Eq. (140) 



in 



^ = 2(+=)iV;(lnS/A). 



The cascading factor here can be presented as 

"^ dE„ ... asjE, 
En "■"' 2tt 






K{'^4 



UNc^^^^^nJIu^ 



A 



(152) 



(153) 
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One can easily see that for the radiative two-jet events the emission pattern is 
given by the qq sample Lorentz-boostcd from the quark cms to the lab system 
(i.e. the cms of qqj), and the corresponding particle multiplicity should be 
equal to that in e+e^ — > qq at W^^ — {pq + pq)^ . 

Now we consider the three-jet event sample when a hard photon is replaced 
by a gluon gi. For a given qqgi configuration the particle flow can be presented 
as 



dQn2 



= — N'{Yq+ , r,) + — N'AY,^ , Yq) 



a+ 



0-1 



'1+ 



h-- l- 



2Cf 



where, in addition to the definitions in Eq. (149) one has 



K{Y), 



(154) 



Y, . mf . 



Y,i 



In 



Eg^/^^ 



A 



(155) 



This formula accounts for both types of coherence: the angular ordering inside 
each of the jet s an d the collective nature of the inter-jet flows. The first three 
terms in Eq. ( |l54| ) are collinear singular as &i ^ and contain the factors 
N' , describing the evolution of each jet initiated by the hard emitters q, q and 
gi. The last term accounts for the interference between these jets. It has no 
collinear singularities an d co ntains the common factor N'(Y,Y) independent 
of the direction n2- Eq. (|l54| ) predicts the energy evolution of particle flows in 
qqg events. 



As follows from ( 148 )-(154) when one is replacing a hard photon by a 
gluon gi with otherwise identical kinematics, an additional particle flow arises: 



/'8TTdN\ 



SirdNgg 

d^n2 



hrdNgg 

d^n2 



(156) 



ai 



KiYguYg) + [h+ + /i_ - I+-]N'{Y). 



Note that for the case of large radiation angles the two cascading factors N' 
become approximately equal and one has, c.f. Eq. (138) 



/87rdiV\ 



(1+) + (1-) - (+=)J KiY)- 
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(157) 



An instructive point is that this expression is not positively definite. One 
clearly observes the net destructive interference in the region between the q 
and q jets. The soft radiation in this direction proves to be less than that in 
the absence of the gluon jet gi. 

7.4 Topology Dependence of Three-Jet-Event Multiplicity 

Now we turn to the connection between the average particle multiplicities in the 
two-jet and three-jet samples of e^e~ annihilation. The particle multiplicity 
in an individual quark jet is formally defined from the process e"*'e~ -^ qq ^ 
hadrons by 



K+.-m = 2N^\E) 



l + O'"^^^) 



W = 2E. (158) 



As we have already discussed, when three or more partons are involved in a 
hard interaction the multiplicity cannot be represented simply as a sum of the 
independent parton pieces, but rather it becomes dependent on the geometry 
of the whole ensemble. 

So, the problem arises of describing the multiplicity in three-jet events, 
Nqqg, in terms of the characteristics of the individual q and g jets. The quantity 
Nqqg should depend on the qqg geometry in a Lorentz-invariant way and should 
have a correct limit when the event is transformed to the two-jet configuration 
by decreasing either the energy of the gluon gi or its emission angle. 



The angular integral of Eq. (148) can be easily checked to reproduce the 



total multiplicity. One can write Nqq^ as 

N,,, = J ^ = N,iY,) + N,iY,) + 2ln^l^N:^{Y). (159) 

Now let us transform this formula to the Lorentz-invariant expression. To do 
this we rewrite 

^^qil) = >^ + lna;+(-), a;+(-) = Eg^q)/E (160) 

and use the expression 

Ng{Y,) = N,iY) + \iix+N'^{Y) + O (^ A^g) • (161) 

Then 

N,,, = 2N,{Y) + In ^±^_^±^ 7V^(y) + o{asN) ^ 2N,{Y;_)[1 + 0{as)] 

(162) 
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with 



Y*_ = Y 



In 



2:+x_a+- 



= In 






Here E* is the quark energy in the cms of qq, i.e. the Lorentz-invariant gener- 
ahzation of a parameter of hardness of the process. 

The muhiphcity Nqqg can be written by analogy as 



iVo, 



dNa 



dOif 



In 



Ql+Ql 
2a+- 



N,{Y,) + N,{Y,) + Ng{Yg) 
2Cf 



(164) 
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where Yg = 1" + lna;i. Proceeding as before, one arrives finaUy at the Lorentz- 
invariant resultE3 



Nggg = [2Ng{Y;_) + N g {Yg^ )] (l+O (^) 



\ TT / 



with 



y; = ha. 



(p+ -pOCp- -pi) ^ , pi± 

2(p+-p_)A2 '^2A' 



(165) 



(166) 



where pix the transverse momentum of gi in the qq cms, c.f. Eq. (146). 

Comparing (159) with (165), we see that replacement of a photon 7 by a 
giuon gi leads to the additional multiplicity 



NgiY*) = Ng,g{W) - Ngg^iW), 



(167) 



which depends not on the giuon energy but on it s transverse momentum, i.e. 
on the hardness of the primary process. -[I Eq. (165) reflects the coherent nature 
of soft emission and has a correct limit when the event is transformed to the 
two-jet configuration. Eq. ( |l65f ) can be rewritten in terms of the observed 
multiplicity in e~^e~ collisionslill as 



Nr,. 



N,+,-{2E*) + lr(pii)7V,+, 



(Pi 



(l+Oias)). (168) 



J Within antenna^attern formalism the observation that a giuon adds multiplicity related 
to its Pi— was made.E3 A similar result was derive d in i tho p tring-length approach to the dipole 



pictureJiSIl more details can be found elsewhere 



168 



169 
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In this formula r{pi±) denotes the ratio of multipUcities in gluon and quark 
jets at cms energy W — pi±. The expression for Nqqg can be presented in 
another forr 



N,<ig 



= [Ng{Y,+) + iV,(ri_) + 27V,(r+_) - Ng{Y+^)] [I + 0{as)], (169) 



where Yij = 



In ( \/{pi ■ Pj)/2A'^] = ln{E*j/A). Such a representation deals 

with multiplicities of two-jet events at appropriat e in variant pair energies 
Wij — 2E*j in the (ij) cms frame. Expression ( |l69| ) has a proper limit, 
2Ng{W/2), when the qqg configuration is forced to a quasi-two-jet one, g{8) + 
qq{8), with a suffi ciently small angle between the quarks. The experimental 
tests of Eq. ( 169 ) can be per forme d wi th the tagged heavy quarks, see also 
Ref. 89. Analogously to (168), Eq. (169) could be presented as 



N,<ia 
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-{2E*) 



1 - 



r{E* 
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Ne+e- (2^: 
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(170) 
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iEi 



N,+,- {2El_) 



Until now we have confined ourselves to the so-called 'unbiased' case and 
have not addressed an issue of the three-jet-selection procedure. As well known, 
different approaches to the three-jet selections employ different definitions of 
the ggg-kinematics. Moreover, most jet definitions give three-jets which are all 
biased. 

To gain some insight on how the jet resolution works, let us first consider 
a sample of two-jet qq-evenis selected in such a way that there are no sub-jets 
with pj^ > p_Lcut (within a fc^-clustering scheme with a resolution parameter 
Pi^cut this means that these are only two primary q- and g-jets). If Nqq{L) is 
the 'unbiased' gq-multiplicity, {Nqq{L) = Ng+^-{L) with L = ln(s/A^)), then 
in events with precisely two jets at resolution j/cut — p\cntl ^ there is a rapid- 
it y pl ateau of length ln(l/ycut), see Fig. O and the multiplicity Nqq{L,Kcut) 
; ^|155i|170||T7ll 



Nqq{L,Kcut) ~ Nqq{Kcut) + {L - K^ut) N'jKcut)] 



(171) 



In 



A2 



76 



The first term corresponds to two cones around the q and q jet directions. This 
also corresponds exactly to an unbiased qq system with cms energy p_Lcut- 
The second term describes a central rapidity plateau of width {L — Kcut) — 
In l/i/cut, in which the limit for gluon emission is given by the constraint Kcut- 
This expressi on for a two-jet event can be generalized for a biased multi-jet 



configuration.liijJtLLilli^ 
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Figure 17: Rapidity plateau in two-jet events. 



The average particle multiplicity in the selected two-jet sample is smaller 
than in an unbiased sample. The modification due to the bias is similar to 
the suppression from a Sudakov form factor. It is formally 0{as), but it also 
contains a factor \n^{s/p\). Thus, it is small for large p±- values but it becomes 
significant for smaller p±. This clearly demonstrates that the multiplicity in 
this restricted case depends on two scales, ^/s and p±cut- The p± of an emitted 
gluon is related to the virtual mass of the radiating parent quark. Therefore, 
the two scales \/s/2 and p±cut represent the energy and virtuality of the quark 
and antiquark initiating the jets. 



Though the leading-logarithmic result in Eq. (171) describes qualitatively 
an impact of the bias, for a quantitative analysis an account of the sublead. 



1 f\^ 
ing effects may be required. These subleading corrections were calculatedliii? 

where the MLLA formula for the biased two-jet-event multiplicity was pre 

sented as 



Nqg{L,Kcut) ~ Nqq{Kcut+Cq) + {L~Kcut-Cq)N'jKcut + Cq); 



(172) 



The results from the ARIADNE Monte Carlo are demonstratedliH^ to be in 
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good agreement with the prediction of Eq. (172). 

Let us return now to the three-jet-event case. As was already discussed, in 
the large- iVc limit the emission of softer gluons from a qqg system corresponds 
to two antennas-dipoles which emit gluons independently. If a gluon jet is 
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resolved with the transverse momentum p±, this imposes a constraint on the 
emission of sub-jets from the two antennas. Thus, the con tribution from each 
antenna is determined by an expression hke Eqs. (171,172). 

An essential issue concerns the definition of the jet re solution scale p± for 
three-jet events. In the Lund dipole approachLL^tLiillliHj p± is defined as 



pIlu 



4b+ ■?!){?- -Pi) 



(173) 



c.f. Eqs. (|14q,|l66|). 

These two p^-definitions agree for soft gluons but deviate for harder glu- 
ons. Wh ile p±hu is always bounded by y/s/2 the scale parameter pi± in 
Eq. ( |l66| ) has no kinematic upper limit in the massless case. For hard glu- 
ons p±hu is of the same order as the parent quark virtuality. It was shown 



172 



that the 0(0^) matrix elements are well described if p±hu is chosen as an 
ordering paramete r in the perturbative branching. Therefore, it looks natural 



to assume 
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that the constraint on further emissions is described by the 
identification p±cut — P±Lu- The multiplicity in a qg-dipo\e with restricted 
p± can be described, analogously to the qq-case, as two forward regions and a 
central plateau. 

We conclude this subsection by focusing on the three-jet configurations 
obtained by iterative clustering until exactly three jets remain, without a spec- 
ified resolution scale, where, hence, the constraint on sub-jet p± is described 
by p_Lcut — P_LLu- This allows us to avoid 'biasing' the gluon jet sample, 
which make s thi s selection procedure convenient for extracting the unbi ased 



multiplicity.li^ This strategy is chosen in the recent DELPHI analysis t!_L£ 
while each event is clustered to precisely three jets 



m^ 



174 



By clustering each eve nt to thr ee jets one gets 
three-jet events (c.f. Eqs. (165-168)) 
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for the multiplicity in 



N„^. 



119 



Nqq{Lqq,Kl^^) + - Ngg{ 



(174) 
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(175) 
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An alternative expression for Nqqg is given ai 



N - 



Nqq{L,Kl^u) + -^NggiKLu). 



(176) 



Note that the consistency between Eqs. (174) and ( |176| ) follows from the fact 
that the total rapidity range in the qg and gg-dipoles, 



L„ 



L„ 



In 



2(P+-Pi) 



In 



2(P--Pi) 



A2 A2 

can be described in two different ways by the equalities 

^qq + J-zqa — J-zqa + 1 = 17 + KLu- 



(177) 



(178) 



Let us emphasize that due to QCD coherence the argument in Ngg is pfj_ in 
(174) and pfj_i^^^ in ( |176| ) and not {2pj_)'^ as one would naively expect. The 
difference between the results of Eqs. (174) and (176) arises from the subleading 
corrections (for example, recoil effects) which are not controllable within this 
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framework. 

Formally the bias is related to order as effects, and while it practically 
does not affect the gluon piece it can be important for the qq-contribution. 
Selecting events with comparatively_large values of p± , where the bias is not 
essential, we arrive at the old resultE2l (c.f. Eq. (165)) 



Nqqg 



Nqq{Lqq) + \Ngg{K) 



[l + 0{as)]. 



(179) 



The effect of the bias becomes more and more important for smaller values of 
P±cut- On the other hand,li2y the bias is negligible if one selects gluon system 
recoiling against two quark jets in the same hemisphere. 



7.5 Experimental Studies of Three- Jet Events 

In this subsection we focus mainly on the new experimental evidence of the 
inter-jet color coherence effe cts from the an alvses of the LEP-1 data, also 
discussed in recent reviews 1*^ 'H^^^ 'H^'341*^ ' ^1^ ' "I The main lesson from the recent 
impressive studies is that now we have (quite successfully) entered the stage 
of quantitative tests of the details of color drag phenomena. Owing to the 
massive LEP-1 statistics, one can perform a very detailed study allowing to 
test even some subtle features of the theoretical expectations. 
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Figure 18: Charged particle flow in tjgg-evejita—j in comparison with analytical 
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Let us mention a few n ew fa cts concerning comparison with the analytical 
QCD predictions. DELPHlliiU have performed the first quantitative verifica- 
tion of the perturbative predictionEj for the ratio R^ 



Rry 



Nqq{qqi) 



(180) 



qqg 



of the particle population densities in the interquark valley in the e^e~ 
and e^e~ — > qq^y events. For a clearer quantitative analysis a comparison was 
performed for the F-shaped symmetric events using the double- vertex method 
for the g-jet tagging. The ratio R^ of the charged particle flows in the qq 
angular interval [35°, 115°] was found to be 

R"^^ = 0.58 ± 0.06. (181) 

This value is in fairly good agreement with the expectation following from 
Eqs. (141) and (142) at Nc — 3, for the same angular interval 



R 



th 



0.65iVg 



0.61. 



(182) 



(Note that thi s rat io is slightly larger than the one predicted for Mercedes- 
type events in (143)). The string/drag effect is quantitatively explained by the 
perturbative prediction and the above ratio does not appear to be affected by 
hadronization effects in an essen tial way. 

Another DELPHI resultli^ concerns the analysis of the threefold-symmetric 
e+e^ ^ qqg events using the double-vertex tagging method. It is shown that 
the string/drag effect is clearly present in these fully symmetric events and 
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it cannot be an artefact due to kinematic selections. The aziniuthal-angle 
dependence of particle density in the event plane is shown in Fig. Hq. Quanti- 
tatively, comparing the minima located at ±[50°, 70°], the particle population 
ratio Rg = Nqg/Nqq in the q — g and q — q valleys is found to be 

Rf"^ = 2.23 ± 0.37. (183) 

This is to be compared with the asymptotic prediction Rg = 2.46 for projected 



ojectea 



rates at central angles, whereas for the above angular interval one finds 

Rf w 2.4 (184) 

in good agreement with the experimental value. The number ( |184D was ob- 
tained from the prediction for the full angular distribution^^ (see subsection 



2.3) which is also displayed in Fig. 18 with the normalization adjusted. The 



calculation takes into account both the inside-jet and inter-jet coherence. As 
can be seen, the relative depth of the two valleys between the jets is well repro- 
duced. Also reasonable is the distribution around the gluon jet direction. On 
the other hand, more particles than expected are found within the quark jets. 
In part this is a consequence of using in the calculation the asymptotic value 
4/9 for l/r, which, as we have already discussed, is approached rather slowly 
from above. Despite these shortcomings, the main effect, the ratio Rg of par- 
ticle densities in the valleys between the two types of jets, has been correctly 
predicted by the analytical calculations. 

Another prediction concerns the beam-energy dependence. At a higher 
energy the angular density in Fig. na gets essentially multiplied by an overall 
factor:E3 the particle density in between the jets rises at a rate comparable to 
the density within the jets. 

An instructive test of the inter-jet coherence can be performed when stu dy- 
ing the particle flow in the direction transverse to the 3-jet event plane.L^ 
Consider the radiation of a soft gluon perpe ndicu lar to the plane of a 3-jet qqg 
event. This radiation can be obtained from ( |l38| ) where the contribution from 
one antenna is simply (ij) = 1 — cosG^ . The ratio Rj^ = dN'i^^ / dN'^'^ of soft 
perpendicular radiation in three and two-jet events is obtained asE^ 



Nc 



Ri_ = -^[2-cosei+-cosei_- — 2-(l-cose+_)]. (185) 



In the large Nc limit we have just the superposition of two dipoles with con- 
tribution ^ 1 — cosOiq. Note the limiting cases R± = 1 for coUinear primary 
gluon emission and the proper gg limit Rj_ = Nc/Cf for the configuration 
with parallel qq (6+- — 0) recoiling against the gluon. 
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Figure 19: Multiplicity witbi&ja 30°-cone perpendicular to three-jet-event plane as a fu nction 
of the event topology.lLL3 The curve represents the perturbative prediction (185). 



1 7S 
DELPHItlEI have presented the first results on the particle yield in the 

transverse direction for the ^-shaped symmetric qqg-events. Fig. p^ shows the 

multiplicity within the cone with the fixed opening angle of 30° perpendicular 

to the event plane. The abscissa is the angle Qi between the low-energy jets 

(the non-leading g-jet and the glu on) . The plotted curve corresponds to the 

perturbative prediction from Eq. ( |l85| ) using Oi+ = 0-| = tt — Oi/2 and 

0i_ — 01 with the normalization left free. The rise of particle flow observed for 

increasing opening angle Gi signals an increase of the "effective color charge" 

corresponding to a transition from a qq to a gg type antenna. The experimental 

data prove to be in a quite good agreement with the perturbative predictions, 

providing a new test of the analytical results independent of the hadronization 

models. 

Of special interest is also the measurement of the yield of particles per- 
pendicular to the production plane in the low momentum range.Ell Within 
the perturbative scenario the very soft particle production should be sensitive 
again only to the color charge topology of the primary emitters. It looks quite 
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challenging to test whether the angular dependence of the particle yield given 
by Eq. (185) holds downwards up to small momenta p ^ Qq. 



If one allows for arbitrary three-jet kinematic configurations, new informa- 
tion can be obtained about the evolution of th e evg nt portrait with the varia tion 

and 



17'; 



DELPH 



lit with the yaria , 



of top logv . see previous subsection. ALEPH, 
OPALLLHIJ have convincingly demonstrated that, in agreement with the QCD 
radiophysics, the event multiplicity in three-jet events depends both on the jet 
energies and on the angles between the jets. These results clearly show the 
predicted toplogical depeiidence of jet properties. 



Recently DELPHI! 



173 



have reported the new results on the total charged 
particle mult ip l icity in F-shaped three-jet events. A crucial point in the DEL- 
PHI analysis p " ''^ is that each event is clustered to precisely three jets. This 



between the 
. Recall that 



allows to determine Nggg as a function of the opening angle, Qi 
low-energy jets and to extract the 'unbiased' gt^-multiplicity, Ng 
in symmetric F-events there are only two scale parameters, ^/s = Mz and 8i, 
so the event-topology dependence can be expressed as a function oiOi only, as- 
suming that the leading j et is not the gluon jet, which is true for most cases. A 
comparison is performedlil^l with analytical resul t s ba sed on Eqs. ( 174 , 176). 
The values of Nqq entering the expressions (174, 176 ) via Eq. (172) are ex- 
tracted from the existing data on the energy dependence of iVg+e- ; using the 



standard para metrizations.l 



Fig 20 



'^1 shows the results of comparison with the predictions of Eqs. ( 174 ) 



and (176). In both cases one finds a good agreemen t with the data. DELPHI 
analysis clearly demonstrates that Eqs. ( 174 , 176) fit the data in the wide 
01-range better than the 'unbiased' Eq. (179). 



7.6 Collective Phenomena in High-p±^ Hadronic Reactions 

Color coherence leads to a rich diversity of interference drag-effects in high-pj^ 
multi-jet events in hadronic collisions (see, for example, Refs. 128,146,155,179, 
161,162) and in 7p and DIS processes (see, for example, Refs. 59,180,181). The 
analysis of collective drag phenomena in these processes is considerably more 
subtle than in e+e~ annihilation due to the presence of colored constituents in 
both initial and final states. 

Recall that during a hard interaction, color is transferred from one par- 
ton to another and the color-connected partons act as color antennas, with 
interference effects taking place in the initial or final states, or between the 
initial and final states. Radiation from the incoming and outgoing partons 
forms jets of hadrons around the direction of these colored emitters. The soft 
gluon radiation pattern accompanying a hard partonic system can be repre- 
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Figure 20: Fits of Eqs. (174) (a) and Eq. (175)(b) lo-the data. The full points depict the 

range of the fit, ^''^ 



sented, to leading order in \INr< as a sum of contributions corresponding to 



the color-connected partons. 
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i[[u_y| jj-^ hadronic high-pj_ reactions by varying 
the experimental conditions (triggers) one may separate the dominant par- 
tonic subprocesses and switch from one subprocess to another. Recall also 
that the length and height of the hadronic "plateau" depend here on different 
parameters: the length is determined by the total energy of the collision, and 
the height and the plateau structure depend on the hardness of the process 
governed, as a rule, by the transverse energy-momentum transfer. Thus, in- 
formation becomes available that is inaccessible in e^e^ annihilation where 
both energy and hardness were given by the value of ^/s. Therefore, the high- 
Pj^ hadronic collisions provide us with a very prospective laboratory for the 
detailed studies of the color-related phenomena. 

As a simple illustrative example l et us cnr isider the topology of events, re- 
sulting f roTTi the high-pj^ gg-scattcring ^1 ^ which has been studied recently 



detailtSl 



Xl ■ X2 S 



(a;i,a;2 -- 1). 



(186) 
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In this case the two crossing processes shown in Fig. pi^(a) and Pl|(b) have ap- 
proximately equal probabilities. However, each of them has its own color topol- 
ogy, and therefore, specific particle flows, as schematically shown in Fig. |2llfc') 
and|2T 




Figure 21: Color antennas for two crossing subprocesses of gg-scattering (a,b) and the draw- 
ings of the corresponding hadronic flows (c,d). 



For the subprocess of Fig. El^(a) the soft particle radiation pattern is 



X <^ (14) + (23) + 



1 



Nc ^\ A 



(187) 



2NcCf 



2 [(12) + (34)] - (14) - (23) - (13) - (24) 



where antenna ij is given by Eq. ( |139| ). In full analogy with the discussions 
of the drag effect in subsection 2.2 one may say that the leading contribution 
(the first term in ( |187| )) has the structure of the sum of two independent 
gg-antennas (14) and (23). This fact also can be mimicked by means of the 
topological picture of the Lund string model. Let us emphasize that in our 
case to each contribution (single antenna) a dynamical distribution corresponds 
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which takes into account the effects of cascade multiplication. Furthermore, 
the perturbative approach permits one to control not only the leading color 
contribution, but also the 0{\/Nc) corrections. 

The complete set of antenna pat terns for various parton scattering sub- 
processes for arbitrary Nc is listed in 
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It has been known for a long timel™lu.££l that an especially bright color 



interference effect arises in the case of large- £'t production of color singlet 
objects, for instance, in V+ jet events (with V = j, W^ or Z). The hadronic 
antenna patterns for such processes are entirely analogous to that in the cele- 
brated string-drag effect in e'^e" — > qqg events. 

Rcccnt lv the first (very impressive) data on W+ jet production from 



D0t£fHi£3 have become available. The color-coherence effects are clearly seen 
and these studies have a very promising future. They may play the same role 
for hadron colliders as the important series of results on in ter-j et studies at 
e~^e~ colliders. The quantitative predictions were presentedli2d for color in- 
terference phenomena in the distribution of soft particles and jets in V+ jet 
production at hadron colliders, in particular the Tevatron ^/s = 1.8 TeV pp 
collider. There are two leading-order processes, qq -^ Vg and qg — > Vq. Each 
has its own distinctive antenna pattern. In principle, the antenna pattern could 
be used as a 'partonometer' to identify the dominant scattering process. 

Omitting the CKM factors for clarity, the matrix elements squared for the 
lowest-order processes (for the case V = W) can be written as 



J2\M\' {q{pi)q{p2) ^ W{p:i)g{pi)) = 



El-^l' (9(Pi).9(P2) ~- W{p^)q{pi)) - 



ghw 



ghw 



1 \ 


t^ + u^ + 2sM^ 


N'c) 


tu 


(188) 


t'^ + s 


2 + 2uM^ 





1 t 

Nc ~ -ts 



where s = (pi +p2)'^,t — {pi — pa)^ and u = (pi — ^4)^. In the soft-gluon 
approximation, the corresponding 2 — > 3 matrix elements are 



J2\M\\qq ^ Wg{g)) - g'^Nc [[U] + [24] - -^ [12] j J2\^\'(l^ ^ ^5) 

(189) 
Y.\M\'{qg ^ Wq{g)) = g'.Nc f [12] + [24] - -1^ [14]) ^|Al|2(q.9 ^ Wq) 



Nh 



with 



_ Pi-Pj _ 1 



zj] ^ -; "^ , = -, m (190) 

86 



o '- 








CO 




Qi 




>^ 


"1 






■^ 1.5 


^ y^ 


Q. 


\ / T 


4—1 


H^v /i 


3 


^^ *y''T 


2 


^4.^ 


1— 


^r=^,0,/ 


% 1 






li:i 


• Data 


g 


— MLLA+LPHD 


■:;::3 


- 


0) 


1 1 1 



71/4 7t/2 371/4 



)3 



Figure 22: Comparison of the jct/W multipliq itjr .ratio from data to the MLLA-LPHD 

predictionli2z 



Note that for these processes, the effect of the soft ghion emission is simply 
to multiply the lowest-order matrix elements squared by an overall factor con- 
sisting of three different antennas, defined according to ( |189| ), one of which 
is suppressed in the large Nc limit. This structure is universal for any elec- 
troweak boson -I- jet production, i.e. V — W^, ^",7,7*. 

It is worthwhile to mention that the study of soft radiation pattern ac- 
companying dijct photoproduction can pro vide a, new tool to distinguish the 
so-called direct and resolved mechanisms J2Mlli2i!l Of special interest here is the 



measurement of the soft particle yield perpendicular to the reaction plancE^ 
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Inter-jet color coherence effects were successfully studied by both thuXIIIEl 
and DO Collaborations in pp collisions at the Fermilab Tevatron Collider Ji2il ^"'^H '^ 
The experimental results look rather promising. Thus, the measurements! °^|^ °'^ 
of the spatial correlations between the softer third jet and the second leading- 
Et jet in pp — > 3 jet +X events has clearly demonstrated the presence of 
initial-to-final interference effects in pp interactions. 

In t he j nost direct way the drag-effect signal is extracted in the recent DO 
analysiaird of W+ jet events by comparing the soft particle angular distribu- 
tions around the colorless W boson and opposing leading-_ET jet in the same 
event. In this study the W boson provides a template against which the soft 
particle pattern around the jet can be compared. Such comparison reduces the 
sensitivity to various biases that may be present in the vicinity of both the W 
boson and the jet. 

Inter-jet coherence effects clearly manifest themselves as an enhancement 
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of soft particle radiation around the tagged jet in the event plane (the plane 
defined by the direction of the W boson and the beam axis) relatively to the 
transverse plane when compared with the particle production around the W 
boson. This is illustrated in Fig. g^ which shows DO results for the }et/W par- 
ticle flow ratio as a function of the azimuthal an gle d together with the pertur- 
bative expectation based on MLLA and LPHD.Ii^ The analytical prediction 
proves to be in good quantitative agreement with the data, thus providing new 
evidence supporting the LPHD picture. 

8 Conclusions 

Perturbative QCD proves to be very successful in its applications to multi- 
particle production in hard processes. Still, the problem of the soft limit of 
the theory and of color confinement is not solved. Therefore, at present stage 
hadroproduction phenomena cannot be derived solely from the perturbation 
theory without additional model-dependent assumptions. The complexity of 
the existing popular models requires Monte-Carlo methods to derive their pre- 
dictions. 

In this review we have concentrated on the question of the extent to which 
the semisoft phenomena in multiparticle production in hard processes reflect 
the properties of the perturbative QCD cascades. 

During the last years the experiments have collected exceedingly rich new 
information on the dynamics of hadronic jets — the footprints of QCD partons. 
New QCD physics results from LEP2, TEVATRON and HERA continue to 
pour out shedding light on various aspects of hadroproduction in multi-jet 
events. The existing data convincingly show that the analytical perturbative 
approach to QCD jet physics is in a remarkably healthy shape. 

The key concept of this approach is the hypothesis of local parton- hadron 
duality, according to which the parton cascade can be evolved down to a low 
virtuality scale of the order of the hadronic masses, where the conversion of 
partons into hadrons occurs. Therefore it is the physics of QCD branching 
which governs the gross features of the jet development. 

In particular, the perturbative universality of jets appearing in different 
hard processes is nicely confirmed. Moreover, the data demonstrate that the 
transition between the perturbative and non-perturbative stages of jet evolu- 
tion is quite smooth, and that the bright color coherence phenomena success- 
fully survive the hadronization stage. 

LEPl proves to be a priceless source of information on the QCD dynam- 
ics. It has benefited from the record statistics and the substantial lack of 
background. We have learned much interesting physics, but there is still a 



need for further refined analyses of the data recorded at LEPl. 

The programs of QCD studies at the LHC and at a future linear e+e" 
collider look quite promising. The semisoft QCD physics steadily becomes 
one of the important topics for investigation in the TEVATRON and HERA 
experiments. 

Concluding this review let us emphasize once more that, of course, there is 
no mystery within the perturbative QCD framework. One is only supposed to 
perform the calculational routine properly. So it is entirely unremarkable that 
the quantum mechanical interference effects should be observed in the pertur- 
bative results. Of real importance is that the experiment demonstrates that 
the transformer between the perturbative and non-perturbative phases acts 
very smoothly. This message could (some day) shed light on the mechanism 
of color confinement. 
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